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ABSTRACT 


The  purpose  of  this  thesis  is  to  develop  a theory  of  the 
mechanical  properties  of  noncross-linked  polymers  at  temperatures 
above  the  glass  transformation.  The  relationship  between  this  theory 
and  the  dielectric  dispersion  of  polar  polymers  is  also  discussed. 

If  the  bond  lengths  and  bond  angles  of  the  molecular  chains 
are  constant  during  extension  of  a polymer,  the  configuration  of  an 
entire  molecule  at  any  time  t may  be  specified  by  the  azimuth  angles 
of  the  bonds  in  the  chain,  and  by  the  two  angular  coordinates  describ- 
ing the  rotation  of  the  chain  around  its  center  of  mass.  The  state 
of  the  material  is  described  statistically  by  f , the  relative  number 
of  molecules  having  each  of  the  possible  configurations. 

For  a noncross-linked  polymer,  the  only  intermolecular 
forces  are  viscous  in  nature,  and  the  force  on  a single  element  (one 
skeletal  atom  together  with  its  side  groups)  is  proportional  to  the 
velocity  of  that  element  relative  to  the  surrounding  atoms.  If  the 
polymer  is  incompressible,  and  if  the  strain  is  everywhere  the  same, 
it  follows  that  the  molecular  configurations  during  extension  of  the 
polymer  may  be  described  by  Smoluchowski's  diffusion  equation: 

V'Pj^.fW/tT]  - |L, 

where  D is  the  diffusion  tensor  of  a molecule,  k is  Boltzmann's  constant, 
and  T is  the  absolute  temperature.  The  potential  V under  which  the 
diffusion  takes  place  is  equal  to  the  rate  of  extension  multiplied  by 


a quadratic  function  of  the  positions  of  the  elements  relative  to  the 
center  of  mass  of  the  molecule. 

Once  D and  the  rate  of  extension  are  known,  the  dif 'usion 
equation  can  be  solved  for  f,  the  probability  of  each  molecular  con- 
figuration, From  f it  is  possible  to  find  the  stress  on  the  polymer. 

In  other  words,  the  diffusion  equation  leads  to  a general  relationship 
between  the  stress  and  the  extension. 

The  solution  of  the  diffusion  equation  includes  the  theoretical 
dependence  of  the  stress  upon  the  temperature.  The  predicted  temperature 
effects  agree  with  the  observed  mechanical  properties  of  high  polymers 
such  as  polyisobutylene. 

If  it  is  assumed  that  the  molecules  in  tho  unplasticized  polymer 
move  in  the  same  way  as  in  a dilute  solution,  the  diffusion  tensor  is 
identical  with  the  one  given  by  Kirkwood  ar.d  Fucss  in  their  treatment 
of  the  dielectric  dispersion  of  polar  polymers  [J.  G.  Kirkwood  and 
R.  H.  Fvoss,  _J.  Chem.  Phys . j),  329  (1941)].  The  stress  relaxation 
calculated  from  this  assumption  decays  too  rapidly  with  time  to  fit 
the  experimental  properties  of  polyisobutylene.  This  discrepancy  is 
believed  to  be  due  to  the  chain  entanglements  which  occur  in  the 
unplasticized  polymer.  However,  the  diffusion  tensor  can  be  corrected 
arbitrarily  so  that  it  gives  the  proper  time  dependence  of  the  stress 
relaxation  and  therefore  includes  the  effect  of  the  chain  entanglements. 
A diffusion  tensor  is  found  which  correctly  describes  the  mechanical 
properties  of  polyisobutylene  over  nine  decades  of  time.  The  corre- 
sponding theoretical  dependence  upon  the  molecular  weight  does  not 
quite  agree  with  experimental  results. 


Ixi  determining  the  response  of  a polar  polymer  tc  an  electric 
field,  Kirkwood  and  Fuoss  use  the  diffusion  tensor  applicable  to  dilute 
solutions.  It  is  possible  to  introduce  into  thair  calculations  the 
new  diffusion  tensor  obtained  from  the  theory  of  extension.  Except 
for  the  molecular  weight  dependence,  the  results  of  this  correction 
agree  reasonably  well  with  the  experimental  dielectric  dispersion  of 
unplasticized  polyvinyl  chloride.  It  therefore  seems  possible  to  find 
a diffusion  tensor  correctly  describing  the  time  dependence  of  a polymer 
response  to  both  mechanical  and  electrical  forces. 
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I.  INTRODUCTION 

An  elastomer  may  be  defined  as  a cross-linked  or  noncross- 
linked  linear  high  polymer  in  which  the  bonds  of  the  molecular  chains 
are  free  to  rotate,  or  at  least  are  free  to  assume  easily  a variety 
of  positions.  Contiguous  molecules  are  also  free  to  move  relative  to 
each  other.  Becauso  the  molecules  lose  much  of  their  freedom  of  motion 
at  temperatures  below  the  brittle  point,  or  glass  transformation,  it  * 
is  necessary  to  specify  the  temperature  range  in  which  a polymer 
behaves  as  an  elastomer.  Examples  of  elastomers  are  polyisobutylene 
at  room  temperature  and  polymerized  sulfur  at  higher  temperatures. 

If  not  worked  excessively,  elastomers  Eire  amorphous,  To  obtain  a 
crystalline  x-ray  pattern  for  polyisobutylene,  for  example,  the  polymer 
must  be  extended  rapidly  at  least  1000)i  [1], 

From  the  molecular  properties  of  these  materials,  several 
things  can  be  inferred  about  their  response  to  mechanical  stress.  When 
no  external  forces  Eire  acting  on  the  polymer  molecules,  they  reach  an 
equilibirum  state  in  which  the  bonds  are  randomly  distributed  among 
their  possible  relative  positions.  On  the  average,  each  molecule  is 
then  coiled  up,  with  the  direct  distance  from  one  end  of  the  molecule 
to  the  other  proportional  to  the  square  root  of  the  molecular  weight  [2], 
External  mechanical  forces  can  extend  a molecule  to  many  times  its 
equilibrium  end-to-end  length  with  only  a slight  change  in  the  internal 

1.  C.  S.  Fuller,  U.  J.  Frosch,  and  N.  R.  Pape,  jJ.,  Am.  Chem.  Soc . 62. 

1905  (1940),  ~ ' 

2.  L.  R.  G.  Treloar,  "Physics  of  Rubber  Elasticity,"  Oxford,  Oxford 

University  Press,  1949,  Chapter  III. 
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energy.  Consequently,  elastomers  are  characterized  by  a small  modulus 
and  a very  large  extensibility. 

The  first  successful  application  of  the  molecular  description 
of  elastomers  was  the  theoretical  treatment  of  the  mechanical  properties 
of  cross-linked  polymers  in  thermodynamic  equilibrium  [3],  If  the 
molecules  are  cross-linked,  as  in  vulcanized  rubber,  the  chains  form 
a network  running  throughout  the  material.  When  stretched,  this  net- 
work can  support  a force  so  long  as  the  cross-links  are  not  disrupted 
chemically.  To  compute  the  force  necessary  to  sustain  a given  strain, 
the  chain  configurations  consistent  with  that  strain  are  counted.  The 
entropy  and  the  force  are  then  computed  by  ordinary  thermodynamic 
arguments.  In  this  calculation,  changes  in  the  internal  energy  are 
neglected,  a reasonable  assumption  in  view  of  tne  fact  that  the  average 
inter -atomic  distance  is  affected  very  little  liy  the  strain.  The 
dependence  of  the  stress  upon  both  the  strain  and  the  temperature, 
calculated  in  this  way,  agrees  reasonably  well  with  experiment  [4], 

Noncross-linked  elastomers,  on  the  other  hand,  cannot  support 
any  stress  in  equilibrium.  The  chains  are  free  to  slip  past  one  another 
and  to  resume  unstressed  configurations.  If  the  material  is  suddenly 

3.  From  a large  literature  one  might  cite: 

K.  H.  Meyer,  G.  von  Susich,  and  E.  Valko,  Kolloid-Z.  59,  208  (1932), 
W.  T.  Busse, _J.  Phys.  Chem.  36,  2862  (1932), 

E.  Karrer,  Phys,  Roy.  39,  857  (1932), 

H,  M.  James,  and  E.  Guth,  J,  Polymer  Sci.  4.,  153  (1949), 

F.  T.  Wall,  J.  Chem.  Phys.  11,  527  (1943), 

P.  J.  Flory,  and  J.  Rehner , Jr.,  J.  Chem.  Phys.  11,  512  (1943), 

W.  Kuhn.  Kolloid-Z.  76,  258  (19367,  and 

L.  R.  G,  Treloar,  Trans . Faraday  Soc . 40T  59  (1944). 

. P.  J,  Flory,  Chem.  Revs.  35,  51  (1944  )j  and  P,  J.  Flory,  N.  Rabjohn, 
and  M.  C.  Shaffer,  <J.  Polymer  Sci.  4 . 225  (1949). 
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extended  and  then  kept  at  a constant  length,  the  stress  decays  very 
slowly  to  zero  from  its  initial  value.  During  a large  part  of  this 
relaxation,  the  force  is  approximately  a linear  function  of  the  log- 
arithm of  the  time  [5], 

The  effect  of  temperature  on  the  stress  relaxation  is  twofold. 
First,  the  force  at  small  times  is  proportional  to  the  absolute  tempera- 
ture, similar  to  the  equilibrium  forces  for  cross-linked  polymers. 
Secondly,  raising  the  temperature  makes  the  relaxation  more  rapid. 
Reducing  the  molecular  weight  also  increases  the  rate  of  relaxation, 
but  does  not  alter  the  magnitude  of  the  force  at  small  times. 

When  an  elastomer  is  polar,  it  not  only  has  interesting 
mechanical  properties  but  it  also  exhibits  an  anomalous  dielectric 
behavior  [6],  The  dielectric  dispersion  of  a polymer  such  as  polyvinyl 
chloride  depends  o:i  the  frequency  of  the  applied  voltage,  and  nas  a 
broad,  low  maximum  in  the  audio  frequency  range,  Debye's  single  relax- 
ation time  does  not  explain  the  frequency  dependence  of  the  loss,  and 
a distribution  of  relaxation  times  is  necessary  to  describe  the  exper- 
imental results.  The  frequency  of  maximum  loss  varies  with  the 
temperature,  indicating  a temperature -time  relationship  similar  to  that 
observed  for  the  mechanical  properties, 

A theoretical  treatment  of  the  dielectric  dispersion  has  been 
given  by  Kirkwood  and  Fuoss  [7],  In  this  theory,  it  is  assumed  that  the 

5.  R.  D.  Andrews,  N.  Hofman-Bang,  and  A.  V.  Tobolsky,  J.  Polymer  Sci. 
J,  669  (1948). 

6,  A summary  of  the  dielectric  relaxation  of  polymers  is  given  b\ 

W.  Kauzmann,  Rev.  Mod.  Phys . 14,  12  (1942). 
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J,  G.  Kirkwood,  and  R,  M.  Fuoss,  J,  Chem.  Phys . j),  329  (1941). 
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the  bond  lengths  and  bond  angles  of  each  polymer  molecule  are  always 
constant,  but  that  the  azimuth  angles  of  the  bends  are  continually 
rotating  due  to  the  thermal  motion  of  the  atoms.  In  the  presence  of 
an  electric  field,  the  azimuth  angles  tend  to  rotate  so  that  the  indi- 
vidual dipoles  of  the  molecule  are  aligned  as  much  as  possible  with  the 
direction  of  the  field.  The  only  forces  opposing  this  change  in  the 
molecular  configuration  are  viscous  drag  forces  on  each  of  the  polymer 
atoms.  The  solution  of  tho  diffusion  equation  describing  this  process 
indicates  that  a molecular  chain  moves  in  sections  of  various  length, 
with  each  section  acting  as  a single  mechanical  and  electrical  unit. 

One  might  say,  roughly,  that  the  distribution  of  possible  section 
lengths  for  a single  molecule  produces  the  distribution  of  relaxation 
times.  This  theory  predicts  a broad  maximum  for  the  dielectric  lo3s 
at  low  frequencies,  but  the  calculated  maximum  is  too  large  and  too 
narrow  to  agree  quantitatively  with  the  experimental  results  for 
unplasticized  polyvinyl  chloride. 

Kirkwood  [8]  has  extended  the  diffusion  theory  to  include  a 
treatment  of  the  nenequilibrium  mechanical  properties  of  cross-linked 
elastomers.  He  assumes  that  the  stress  on  the  material  is  carried 
entirely  by  the  cross-links  between  the  polymer  chains.  Each  molecule 
tends  to  increase  in  length  under  the  stress -dependent  force  across 
its  tiepoints,  with  the  motion  of  the  molecule  being  retarded  by  the 
viscous  forces  on  each  of  the  polymer  atoms.  The  results  of  this  theory 
have  not  been  compeared  with  the  experimental  properties  of  elastomers. 


8,  J.  G.  Kirkwood,  _J.  Chem.  Phvs.  14,  51  (1946), 


5 


In  Kirkwood's  theory,  it  is  assumed  that  the  viscous  force 
on  an  atom  in  a given  molecule  is  proportional  to  the  velocity  of  that 
atom  relative  to  the  center  of  mass  of  the  molecule.  Actually,  however, 
the  viscous  force  depends  on  the  motion  of  the  atom  relative  to  its 
immediate  surroundings.  That  is  to  say,  the  force  is  proportional  to 
the  velocity  of  the  atom  minus  the  velocity  of  the  material  around  the 
atom,  where  both  motions  may  be  measured  relative  to  the  center  of  mass 
of  the  molecule.  The  velocity  of  the  material  at  any  point  is  zero 
during  the  application  of  an  electric  field,  but  is  not  zero  if  the 
polymer  body  is  being  strained  by  a mechanical  force.  For  extension, 
it  may  be  shown  readily  that  this  velocity  is  proportional  to  the  rate 
of  extension  of  the  polymer  body.  Consequently,  there  is  a term  in 
the  viscous  force  which  is  not  related  directly  to  the  stress  on  the 
polymer  or  to  the  motion  of  the  molecule,  but  depends  instead  on  the 
rate  of  extension.  This  term  is  neglected  by  Kirkwood  although  it 
can  contribute  to  the  intermolecular  forces. 

The  importance  of  this  term  may  be  seen  most  clearly  if  we 
consider  a noncross-linked  polymer . In  Kirkwood's  theory,  an  external 
stress  cannot  be  applied  to  a polymer  molecule  unless  it  is  cross- 
linked  to  the  rest  of  the  material.  This  comes  about  because  it  is 
assumed  that  the  stress  is  carried  entirely  by  the  cross-links.  The 
predicted  modulus  of  a noncross-linked  elastomer  is  therefore  zero,  a 
result  that  does  not  agree  with  experiment. 

In  this  thesis,  the  mechanical  properties  of  an  elastomer 
will  be  derived  by  a method  which  includes  the  effects  of  the  viscous 
intermolecular  forces..  We  shall  begin  by  relating  the  rate  of  extension 
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of  the  polymer  to  the  flow  of  the  material  around  the  center  of  mass 
of  a given  molecule.  The  way  in  which  the  flow  affects  the  viscous 
forces  or.  tho  molecule  will  then  be  derived.  A diffusion  equation  will 
be  found  which  describes  +he  way  in  which  the  molecules  change  their 
configurations  under  these  forces.  Only  noncross-linked  polymers  will 
be  considered  in  order  to  emphasize  the  strain-dependent  nature  of  the 
ir.termolecular  forces. 

Qualitatively,  the  diffusion  equation  indicates  that  the 
molecules  move  in  the  following  way.  If  there  are  no  external  forces 
on  the  polymer,  the  thermal  motion  of  the  molecules  tends  to  orient 
them  randomly,  and  the  distribution  of  bond  directions  is  spherically 
symmetric.  During  extension  of  the  polymer,  the  azimuth  angles  rotate 
so  that  the  molecules  tend  to  line  up  in  the  direction  of  the  extension. 
If  the  elastomer  is  then  kept  at  a constant  length,  the  molecules  slowly 
return  by  their  thermal  motion  to  a spherically  symmetric  distribution. 
This  latter  process  is  accompanied  by  the  stress  relaxation. 

Once  the  conf ugurations  of  the  molecules  arc  known,  one  can 
calculate  the  stress  on  the  elastomer.  Since  the  configurations  are 
determined  by  the  rate  of  strain,  it  is  possible  to  find  a general 
relationship  between  the  stress  on  the  elastomer  and  its  resultant 
extension.  This  relationship  will  be  expressed  in  terms  cf  a distri- 
bution of  mechanical  relaxation  times. 

The  theoretical  stress  relaxation  obtained  in  the  above  manner 
has  the  correct  temperature  dependence,  but  decays  too  rapidly  in  time 
to  agree  with  the  experimental  data  for  unplasticized  polyisobutylene. 
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It  is  believed  that  the  primary  reason  for  this  discrepancy  is  that  the 
motion  of  a polymer  molecule  is  altered,  by  its  entanglements  with  the 
neighboring  chains  These  entanglements  retard  the  changes  in  the  con- 
figurations of  the  molecule  and  effectively  increase  the  resistance 
constant  of  each  atom,  where  the  resistance  constant  is  the  ratio  of 
the  viscous  force  on  an  atom  to  its  velocity.  In  other  words,  the 
viscous  force  on  an  atom  in  an  unplasticized  polymer  will  be  much  greater 
than  the  corresponding  force  in  a dilute  polymer  solution.  Furthermore, 
the  magnitude  of  the  force  will  depend  upon  the  position  of  the  atom 
in  its  polymer  chain.  In  order  to  introduce  the  effect  of  the  entangle- 
ments into  the  theory,  the  resistance  constants  of  the  atoms  will  be 
made  an  arbitrary  function  of  distance  from  the  ends  of  the  chain, 
contrasting  with  Kirkwood's  assumption  that  the  resistance  constants 
are  tho  same  for  all  the  atoms.  A set  of  constants  will  be  found  which, 
when  substituted  into  the  diffusion  equation,  gives  the  correct  depend- 
ence of  the  stress  relaxation  upon  the  time  as  well  as  upon  the  temperature. 
The  dependence  upon  the  molecular  weight  is  not  predicted  correctly. 

The  treatment  of  the  dielectric  dispersion  oy  Kirkwood  and 
Fuoss  is  applicable  only  to  dilute  solutions  of  a polar  polymer,  because 
it  does  not  include  any  interactions  between  the  polymer  molecules. 

However,  the  theory  can  be  corrected  for  the  chain  entanglements  by 
introducing  the  set  of  resistance  constants  which  gives  the  correct 
mechanical  properties  of  an  unplasticized  polymer.  This  correction 
will  be  carried  out  in  the  last  section  of  this  thesis.  The  resultant 
dielectric  dispersion  agrees  quantitatively  with  the  experimental  prop- 
erties of  unplasticized  polyvinyl  chloride. 
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Except  for  the  molecular  weight  dependence^  it  would  therefore 
appear  that  the  ad  hoc  set  of  resistance  constants  correctly  describes 
the  molecular  diffusion  within  an  unplasilcized  polymer  at  temperatures 
above  its  glass  transformation. 
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II.  IKE  MOLiJjULAR  PROCESSES  RESULTING  FROM  MECHANICAL  EXTENSION 

This  section  describes  the  effects  of  mechanical  extension 
upon  the  positions  and  configurations  of  the  molecules  in  a noncross- 
linked  elastomer.  By  evaluating  the  intermolecular  forces  which  result 
from  extension,  we  can  find  the  diffusion  equation  satisfied  by  the 
molecular  configurations.  It  is  then  shown  that  solution  of  the  dif- 
fusion equation  leads  to  a relationship  between  the  stress  on  the  polymer 
and  the  resultant  strain.  This  relationship  includes  the  dependence 
of  the  stress  and  the  strain  upon  the  time. 

A.  The  Mac-rosocpic  Flow  During  Extension 

Consider  a right  circular  cylinder  of  elastomer,  of  unit 
cross-section  and  unit  length  when  unstressed.  The  dimensions  of  the 
cylinder  are  much  larger  than  the  dimensions  of  the  molecules.  Let 
the  cylinder  be  extended  parallel  to  its  axis  by  a stress  s(.t)  applied 
to  the  ends  of  the  cylinder.  The  stress  is  an  arbitrary  function  of 
the  time  t . The  macroscopic  strain  e(t)  is  the  relative  increase 
in  the  length  of  the  cylinder.  It  will  be  assumed  that  £(t)  is  always 
much  less  than  unity  and  that  the  area  of  the  cylinder  changes  only 
slightly  during  extension.  It  will  also  be  assumed  that  any  heat  gen- 
erated by  the  application  of  the  external  force  to  the  body  is  carried 
off  rapidly  enough  so  that  the  temperature  T of  the  elastomer  is  constant. 

Let  us  erect  a coordinate  system  at  the  center  of  mass  of  any 
one  of  the  molecules  in  the  cylinder,  with  the  h axis  parallel  to  the 
axis  of  the  cylinder.  During  extension  of  the  polymer,  the  origin  of 
this  coordinate  system  moves  with  the  center  of  mass  of  the  molecule. 

! 
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The  usual  cylindrical  coordinates  p , , ai id  * will  be  used  to  measure 

position  relative  to  the  coordinate  system. 

Consider  a volume  which  has  dimensions  small  compared 

with  the  dimensions  of  the  elastic  body  under  consideration,  but  large 
compared  with  tne  size  of  a single  atom.  Also,  let  us  define  a time 
interval  A"t  which  is  small  compared  with  the  times  involved  in  any 
experiment  on  the  polymer  body,  but  large  compared  with  the  time  inter- 
val-- necessary  to  describe  any  single  atomic  collision  or  interaction. 
Let  v denote  the  average  velocity  of  all  the  atoms  inside  AV  during 
the  time  interval  At  . This  velocity  is  measured  relative  to  the 
coordinaie  system  moving  with  the  molecule.  It  will  be  assumed  that 
AV"  and  At  exist  such  that  Y is  a continuous  and  differentiable 
function  of  time  and  of  position  in  the  body. 

If  the  strain  inside  the  cylinder  is  everywhere  the  same, 
the  velocity  v is  independent  of  the  choice  of  molecule  upon  which 
the  coordinate  system  is  erected.  From  this  fact,  an  expression  may 
be  found  for  y an  terms  of  the  rate  of  extension.  By  cylindrical 
symmetry,  v will  not  have  a component  in  the  y direction,  nor  will 
any  of  its  components  be  a function  of  ^ , Also,  it  may  be  shown 
readily  that  the  a component  of  -v-  is  everywhere  independent  of  p , 

and  that  the  p component  of  \r  is  independent  of  s.  . We  can  summarize 

ft 

these  results  by  the  vector  equation 

V3*  v ~ °- 

it 


V3  denotes  the  usual  three-dimensional  vector  operator  (see  page  16). 


1 


i 

i 

j 
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A second  property  of  v*  is  obtained  from  conservation  of  the 
mass  of  tile  body.  Tf  the  elastomer  is  incompressible,  the  density  is 
independent  of  the  time  and  the  divergence  of  nr  *s  zero: 

V • v-  = 0.  (2) 


Finally,  the  boundary  conditions  upon  v*  may  be  given.  If 
e(t)  is  much  less  than  unity,  the  limit  of  (' as  r becomes  very 
large  is 


f Y)2  “ 2 dt  ‘ (3) 

X °o 

At  the  origin  of  the  coordinate  system,  \r  is  identically  zero. 

Equations  (l)  and  (2),  together  with  these  boundary  conditions,  uniquely 
define  the  velocity  function.  The  r component  of  the  velocity  is  given 
everywhere  by 


(v)  » Z 


dc 

dt 


(4) 


and  the  p component  by 

= ~ J P dt  ‘ (5) 

From  these  equations,  we  may  show  readily  that  the  trajectories 
of  the  material  relative  to  the  molecular  coordinate  system  are 

p z = constant.  (g  ) 

This  set  of  trajectories  implies  that  every  small,  right  circular 
cylinder  of  material,  with  its  axis  parallel  to  the  z axis,  always 
moves  inside  the  main  cylinder  so  as  to  retain  its  cylindrical  shape 
and  constant  volume.  The  flow  pattern  of  the  material  is  shown  in 

Fi  iriir  n 1 
* -o  — ^ - • 
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B„  Description  of  the  Molecular  State 

To  define  the  state  of  the  polymer  at  any  time  t , it  is 
necessary  to  add  to  y a description  of  the  configurations  of  the  mole- 
cules. The  instantaneous  configuration  of  a molecule  relative  to  a 
coordinate  system  erected  at  its  center  of  mass  can  be  specified  by  the 
direction  of  the  center  bond  and  by  the  azimuth  angles  along  the  polymer 
chain.  If  de/it  is  small,  the  bond  lengths  and  bond  angles  arc  constant 
during  extension.  The  side  groups  of  the  molecule  may  be  free  to  change 
position  internally  and  to  rotate  around  the  bonds  connecting  them  to 
the  main  chain,  but  these  motions  can  be  neglected  if  the  side  chains 
are  relatively  short. 

The  zenith  and  azimuth  angles  © and  d will  be  used  to  specify 
the  direction  of  the  center  bond.  If  the  molecular  chain  consists  of 
2r\-l  single  bonds,  the  azimuth  angles  can  be  denoted  by  X3,X'3  , 

numbered  from  one  end  of  the  chain.  Exact  definitions  of  B J and 
^2 ^an-<  are  Siven  *n  Appendix  I.  The  configuration  of  an  entire 
molecule  is  specified  by  its  X2  , X2n_(  . 

Without  loss  in  generality,  it  may  be  assumed  that  all  of  the 
molecules  have  tha  same  structure  and  molecular  weight,  and  are  indis- 
tinguishable. Consequently,  the  state  of  the  molecules  at  any  time  ~t 
may  be  specified  by  the  number  of  molecules  having  each  possible  config- 
uration, without  specifying  which  molecules  have  which  configurations. 

For  this  purpose,  ^ ( © , $■.  X2 , . . . , , ) is  defined  such  that 

f sined0d^dXz...dX2h_(  is  the  fraction  of  the  molecules  having 
configurations  between  © , <f> , X2  , X^  t .. . , ( and  0 + d©,  4>+d<f>, 

d , . . . , ( + d X2  . The  number  of  molecules  is  constant. 
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and  the  integral  of  f over  all  possible  configurations  is  always  unity: 

\S.6%=  '•  (7) 

The  configurations  of  the  molecules  may  be  described  most 
conveniontly  in  the  Zt\  dimensional  coordinate  space  in  which 
©,  •>  X-2  , • • • , X-2n  _ | are  orthogonal  coordinates.  The  differential 

volume  element  dc^  in  this  space  is  equal  to  sin©d©el^  . : . <d  • 

If  the  configurations  of  all  the  molecules  of  the  polymer  body  are  plotted 
as  points  in  this  space,  the  configurational  probability  ^ i3  propor- 
tional to  the  density  of  the  points. 

The  velocity  of  all  the  points  in  the  volume  element 
when  averaged  over  the  time  interval  Ht  , is  denoted  by  o . 
components  of  u are  9,  $ sin©,  Xz  , , X-4  , . . . , X2n|  . The 

velocity  u at  a particular  point  in  the  space  gives  the  rate  of  change 
of  the  corresponding  molecular  configuration.  Specification  of  u over 
the  entire  space  describes  the  way  in  which  all  of  the  configurations 
are  changing  with  time. 

In  summary,  the  state  of  an  incompressible  elastomer  may  be 
specified  for  our  purpose  by  £ , the  probability  of  each  configuration. 
The  rate  at  which  the  state  is  changing  may  be  specified  by  the  velocities 
y and  o . By  using  this  statistical  description  of  the  molecular 
state,  we  can  find  the  response  of  an  elastomer  tc.  mechanical  and 
electrical  forces, 

C,  The  Diffusion  Equation  for  the  i-iolecular  Configurations 

i 


Once  the  state  of  the  elastomer  is  defined,  it  is  possible 
to  describe  the  inter molecular  forces.  Knowing  these  forces,  we  can 
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obtain  the  equation  of  notion  of  the  molecules.  Because  the  stress  and 
strain  are  measured  as  averages  over  times  at  least  as  largo  as  the 
interval  /it  , the  forces  acting  on  a given  molecule  may  also  be  aver- 
aged over  At  . Only  that  average  force  need  be  considered  in  finding 
an  equation  of  motion. 

If  the  side  groups  are  small,  the  average  motion  of  the  groups 
on  any  single  atom  of  the  main  molecular  chain  is  the  same  as  the  average 
motion  of  the  atom  itself.  Hence,  the  entire  molecule  including  its 
side  groups  may  be  divided  into  2 o equal  elements  which  may  be  treated 
as  rigid  bodies.  Each  element  is  at  a constant  bond  length  and  bond 
angle  from  the  next  element  along  the  chain,  but  is  free  to  rotate  around 
the  azimuth  angle  of  the  bond.  The  forces  on  each  element  consist  of 
the  constraint  forces  required  to  keep  the  bond  lengths  and  bond  angles 
constant,  plus  the  average  external  force  due  to  the  surrounding  atoms. 
The  constraint  forces  need  not  be  considered  if  the  molecular  motion 
is  described  in  the  internal  coordinate  space  of  6,  <f> , • 

The  thermal  motion  of  an  element  consists  of  jumps  from  one 
stable  position  to  another,  accompanied  by  rearrangement  of  the  surround- 
ing atoms.  Let  y ; denote  the  velocity  of  the  i-th  element  in  the 
chain  (numbered  from  one  end).  This  velocity  is  measured  relative  to 
the  center  of  mass  of  the  molecule,  and  is  averaged  over  the  time 
interval  At  . If  the  energy  of  thermal  agitation  is  large  compared 
with  the  energy  required  for  a single  jump,  many  such  motions  occur 
during  At  , To  a good  approximation,  -y  is  a continuous  function 
of  time.  This  approximation  implies  that  u is  also  continuous. 
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From  this  description  of  the  motion  of  the  molecules,  it  may 
be  postulated  immediately  that  the  average  intermolecular  force  F-  on 
the  i-th  element  is  viscous  in  nature,  and  is  proportional  to  the  velocity 
of  the  element  relative  to  its  surroundings: 


£ = SitY-Yi*.  <8) 

where  v is  measured  at  the  position  of  the  element  in  question.  The 
resistance  constant  ^ is  assumed  to  be  independent  of  the  velocity 
and  position  of  the  element. 

One  part  of  F-  may  be  derived  from  a potential  V-  , defined 


by 


(9) 


As  in  Equation  (8),  y is  measured  at  the  position  of  the  i-th  element. 
From  Equations  (4;  and  (5)  we  may  show  readily  that 


v-  =?  5T  2 ■ <10) 

Introduction  of  V;  into  Equation  (8)  gives 

fi  = <n> 

For  an  unplasticized  polymer,  there  is  no  need  to  consider 
the  hydrodynamic  interactions  between  elements  of  the  same  molecule. 

The  flow  represented  by  y describes  the  motion  of  all  the  atoms  near 
the  i-th  element,  including  those  on  the  same  chain.  therefore 

includes  the  forces  between  elements  of  the  same  molecule.  This  case 
is  quite  different  from  that  of  a dilute  solution  of  a polymer,  in  which 
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the  hydrodynamic  interactions  must  be  considered  separately  from  the 
viscous  forces  due  to  the  flow  of  the  solvent  [9], 

All  of  the  viscous  forces  acting  on  a molecule  are  equivalent 
to  a generalized  force  in  the  internal  coordinate  space  of 

Qt  i *3  > ‘ This  force  is  given  by 

= — VVf  — /o-  U,  (12) 

where  £ is  the  resistance  tensor.  The  components  of  p may  be  found 
from  the  resistance  constants  by  computing  71  , the  torque  around  an 
angle  X produced  by  rotation  around  some  other  angle  X.  : 

T-  03) 

j 

The  calculations  for  p are  given  in  Appendix  I.  The  potential  V ^ is 
the  sum  of  all  the  individual  potentials  of  the  elements: 

yT  u*> 

It  should  be  noted  that  the  V in  Equation  (12)  (and  in  all  succeeding 
equations)  is  the  vector  operator  associated  with  the  internal  coordinate 
space,  and  consequently  has  2 n dimensions.  The  V3  in  Equations  (9) 
and  (ll)  is  the  usual  three-dimensional  quantity. 

If  the  polymer  body  is  in  equilibrium  and  there  are  no  external 
forces  on  the  molecules,  the  angles  Q , r , ~ a > * - • ■>  ^ 2„_,  are  free  to 
assume  any  possible  set  of  values  (assuming  that  the  bond  rotations  are 

3.  F.  Bueche,  J,  Chem.  Phys . 20,  1959  (1952).  The  hydrodynamic  inter- 
actions and  their  effects  in  polymer  solutions  are  described  by 
J.  G.  Kirkwood,  and  J.  Riseman,  _J_.  Chem.  Phys . 16,  565  (1948); 

J.  G.  Kirkwood,  Rec.  Trav.  Chim.  68,  649  (1949);  and 

J.  G.  Kirkwood,  and  P.  L.  Auer,  _J,  Chem.  Phys . 19,  281  (1951). 
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entirely  unrestricted).  Owing  to  the  thermal  motion  of  the  molecules, 
all  poosvble  configurations  of  the  molecules  axe  attained  during  a long 
interval  of  time.  The  configurational  probability  f is  equal  to  £0 , 
a constant  in  the  internal  coordinate  space: 

* = r ^r'.  (15) 

Furthermore,  u is  identically  zero. 

If  the  body  is  not  in  equilibrium,  .f  is  not  equal  to  .f „ . 

The  thermal  motion  of  the  molecules  tends  to  shift  £ gradually  back 
toward  the  equilibrium  distribution,  and  u is  not  zero.  The  effective 
force  (averaged  over  At  ) that  is  equivalent  to  this  tendency  to  return 
to  the  equilibrium  state  is  -!v„T  7f/{  , where  k0  is  Boltzmann's 

constant.  During  this  motion,  ^ must  satisfy  the  equation  of  continuity. 

*•(«)— IS-  «*> 

When  there  are  no  electric  torques  acting  on  the  molecules, 
the  sum  of  the  thermodynamic  force  and  the  viscous  force  is  equal  to 
zero  (assuming  negligible  inertia  for  the  molecules): 

Ft-  k.TVf/(  = O.  (17) 

Introduction  of  the  diffusion  tensor  D , defined  as 

D = )<©TyP  \ (18) 

i 

and  use  of  Equations  (12),  (16),  and  (17)  lead  finally  to  the  diffusion 
equat ion : 

i V-p-[Vf  + |W,/1.T]  = If.  (19) 

i 

i 


I 
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This  equation  can  be  solved  for  £ if  the  potential  and  the  dif- 

fusion  tensor  0 are  known  . 

The  diffusion  equation  may  be  written  in  the  form  of  Equa- 
tion (19)  only  when  exists,  that  is,  when  the  viscous  forces  are 
"conservative".  The  existence  of  a potential  follows  immediately  from 
the  fact  that  the  curl  of  the  velocity  is  zero  (Equation  1).  If  the 
external  stresses  applied  to  the  elastomer  were  shear  forces,  then  the 
curl  of  the  velocity  would  not  be  zero  and  a scalar  potential  would  not 
exist.  The  solution  of  the  diffusion  equation  for  shear  forces  will 

ft 

not  be  studied  in  this  thesis  . 

If  an  electric  field  instead  of  a mechanical  stress  is  applied 
to  the  polymer,  the  only  change  in  the  diffusion  equation  is  the  replace- 
ment of  by  Ve  , the  electric  potential  of  a molecule  in  the  field. 

The  same  diffusion  tensor  is  applicable  in  both  cases.  The  diffusion 
equation  is  then  identical  with  that  used  by  Kirkwood  and  Fuoss  to 
compute  the  dielectric  response  of  a polar  polymer  [7],  In  general, 
for  tlr*  application  of  either  an  electric  field,  a mechanical  extension, 
or  both  s imult aneous ly , the  equation  of  motion  for  ^ may  be  written  as 

v-P'[vf  + fW/y]  - |f-,  (20) 


*This  section  is  not  meant  to  be  a demonstration  of  the  applicability  of 
Gmoluchowski’s  diffusion  equation  to  elastomers  at  temperatures  above 
the  glass  transformation.  It  is  intended,  rather,  to  show  that  the 
applicability  of  the  diffusion  equation  is  a reasonable  assumption,  and 
that  the  potential  for  noncross -linked  molecules  is  then  given  by 
Equation  (14).  For  a discussion  of  Smoluchowski's  equation,  see 
S.  Chandrasekhar,  Rev.  Mod.  Phys . la . 1 (1943). 

rFor  a discussion  of  shear  forces  in  solutions,  see  [9],  A somewhat 
different  method  of  approach  is  given  by  P.  Debye  and  A.  M.  Bueche, 
jJ.  Cham.  Phys,  16,  573  (1948). 
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where  the  total  potential  is 

V = Ve  4 Vf  . (21) 

This  potential  applies  only  to  noncross-lir.Ucd  polymers  and  dees  not 
include  the  forces  on  a molecule  due  to  permanent  cross-links.  No  study 
of  such  forces  will  be  made  here. 


D . The  Relationship  between  the  St  rest',  and  the  Strain 

Retaining  to  the  problem  of  mechanical  extension  of  a noncross- 
linked  elastomer,  we  can  derive  a general  expression  for  the  mechanical 
energy  required  for  extension.  The  rate  at  which  energy  is  withdrawn 
from  the  surrounding  atoms  by  each  molecule  is  the  sum  of  the  scalar 

r~  * ^ 

products  of  r-  and  v j 


(~  ) 

'5t  'molecule 


(22) 


Replacing  -y  by  (v-  y;)  •+•  v-  and  using  Equation  (8), we  have 


(23) 


The  corresponding  equation  in  the  internal  coordinate  space  is 

*«■  = 

From  Equations  (12),  (17),  and  (18),  we  obtain 


(*§) 

’ St  'molecule  l 


(24) 


( — ) 

'St  molecule 


= _ V£_p_w/' 

* 


(25) 


1 y*,  the  complex  conjugate  of  y,  is  introduced  to  ensure  that  the 
energy  is  real. 
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If  the  total  number  of  molecules  per  unit  volume  is  n0  , the 
number  of  molecules  with  configurations  between  X-2  > ■ • • » 

and  6+cjS,  ...  , + is  dcj^  i 

where  dc^  is  $m©dS  <J  ? dX2...  dX2o_(  . The  total  rate  at  which 
energy  is  introduced  into  the  body  is  the  rate  at  which  energy  is  supplied 
to  each  molecule,  multiplied  by  the  number  of  molecules  having  the  same 
configuration,  and  integrated  over  the  internal  coordinate  space: 

H = -D  W(-dV  (26) 

Applying  the  divergence  theorem  and  noting  that  the  coordinate  space 
has  no  boundaries  (it  is  everywhere  re-entrant  as  the  angles 
0,^,  h~i  change  by  TT  or  2n  ),  we  obtain 

H - {w) 

Thus  the  total  rate  at  which  energy  is  supplied  to  the  elastomer  may  be 
calculated  from  the  potential  and  the  configurational  probability.  It 

is  not  necessary  to  determine  explicitly  any  of  the  forces  or  velocities. 
In  order  to  interpret  the  rate  of  energy  transfer  in  physical 
terms.  Equation  (26)  may  be  considered  in  a slightly  different  form: 

= (28) 

* 

- *ATJ  % <29) 


or 


St 


ae 

3t 


I 
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By  substituting  for  V-  D • ^7V£/k0T)  from  the  diffusion  equation 

and  applying  the  divergence  theorem,  it  is  possible  to  obtain: 

| * "-k-T  I jf ' 'A<  dt  + ,30) 


If  the  entropy  per  unit  volume  were  defined  as 

5=  (31) 

the  first  term  in  Equation  (30)  would  be  the  time  derivative  of  TS  . 

It  represents  the  rate  at  which  energy  is  stored  in  the  elastomer , and 
is  dependent  only  upon  the  configurations  of  the  molecules.  Over  a 
complete  cycle  of  sinusoidal  motion,  this  term  does  not  contribute  to 
the  total  energy  loss.  The  second  term  in  Equation  (30)  is  a positive 
definite  quantity  that  increases  the  energy  supplied  to  the  elastomer 
monotonically  with  the  time.  It  represents  the  viscous  loss  due  to 
the  relative  motions  of  the  molecules,  and  is  the  rate  at  v?hich  energy 
is  dissipated  as  heat. 

It  is  now  possible  to  find  an  expression  for  the  stress  s(t). 


The  rate  at  which  energy  is  supplied  to  the  body  is  equal  to 


(32) 


Equating  this  energy  to  the  energy  given  in  Equation  (27),  we  obtain 
finally 


(33) 


This  equation  completes  a general,  description  of  the  phenomena 
occurring  during  extension  of  a noncross-linked  polymer.  If  the  exten- 
sion c(i)  is  given.  Hie  veloci+y  v describing  the  relative  motions  of 
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the  molecules  may  be  calculated  from  Equations  (4)  and  (5).  Equation  (14) 
gives  the  potential  in  terms  of  dc/dt  , If  the  resistance  tensor 

and  its  inverse  are  known,  the  diffusion  equation  may  then  be  solved 
to  find  the  internal  configurations  of  the  molecules.  Substitution  of 
and  £ into  Equation  (33)  gives  the  stress.  In  other  words,  it 
is  possible  to  find  the  external  stress  that  must  be  impressed  on  the 
cylindrical  polymer  body  in  order  to  produce  a given  extension. 

The  above  set  of  equations  is  equivalent  to  a distribution 
of  relaxation  times  defining  the  general  relationship  between  s(t) 
and  e(t)  . In  Section  IV,  a method  will  be  given  for  calculating  the 
distribution  of  relaxation  times  once  the  resistance  tensor  is  known. 


III.  QUALITATIVE  DESCRIPTION  OF  THE  MOLECULAR  MOTION  DURING  EXTENSION 


Before  deriving  a general  solution  for  s(t)  , it  might  be 
pertinent  to  discuss  qualitatively  the  molecular  motion  resulting  from 
various  types  of  extension.  If  the  material  is  initially  in  equilibrium, 
the  molecules  are  randomly  oriented  and  have  equal  probability  of 
assuming  all  possible  configurations.  The  molecules  move  from  one 
configuration  to  another  by  random  thermal  motion,  and  there  is  no  net 
gain  or  loss  in  the  number  of  molecules  having  each  configuration. 

This  state  corresponds  to  ^ = £0  and  u = O . 

Let  us  suppose  that  the  material  is  now  extended  at  a constant 
rate.  The  molecules  move  past  one  ancih-n'  in  a continuous  stream. 
Relative  to  the  center  of  mass  of  any  given  molecule,  the  flow  may  be 
represented  by  the  lines  in  Figure  1.  The  molecules  move  in  from  the 
sides  and  out  toward  the  top  and  bottom  of  the  figure.  This  motion 
tends  to  pull  the  randomly  oriented  molecules  into  alignment  with  the 
direction  of  extension.  Each  element  along  a molecule  is  subjected 
to  a viscous  force  proportional  to  its  velocity  relative  to  the  sur- 
rounding material.  At  first,  when  the  extension  is  small,  the  elements 
are  free  to  move  nearly  at  the  same  velocity  as  the  surrounding  atoms, 
and  the  forces  on  the  molecule  are  small.  As  the  flow  continues,  the 
constant  bond  lengths  and  bond  angles  prevent  the  elements  from  following 
the  flow  lines,  and  the  forces  on  the  molecule  increase.  After  very 
long  times,  a steady  state  of  viscous  flow  is  reached,  in  which  the 
molecules  on  the  average  are  not  changing  their  configurations  although 
they  are  still  moving  relative  to  each  other.  The  external  stress,  an 


?A 


average  of  the  molecular  forces  through  a cross-section  of  the  material, 
has  then  risen  to  a constant  maximum  value.  The  probability  of  each 
configuration  is  again  independent  of  the  time,  but  is  now  given  by 

*fVVk.T=o,  (34) 

- A -V*/K>T 

or  sj.  — A e > (35) 

where  A is  a normalizing  constant.  In  other  words,  the  configurations 
of  the  molecules  are  described  by  a Boltzmann  distribution,  and  are  in 
thermodynamic  equilibrium  even  though  the  macroscopic  dimensions  of  the 
elastomer  are  changing  with  time. 

Let  us  assume  that  the  extension  is  now  suddenly  stopped,  and 
that  the  length  of  the  cylindrical  body  is  fixed.  The  flow  of  the 
molecules  past  one  another  is  zero.  Each  molecular  chain  is  free  to 
diffuse  back  to  the  state  ^ , assuming,  of  course,  that  no  crystal- 

lization has  occurred.  The  contracting  molecules  tend  to  drag  the 
surrounding  atoms  back  with  them  along  the  flow  lines,  and  an  external 
force  is  necessary  to  keep  the  material  at  constant  extension.  After 
a very  long  time,  f again  returns  to  f0  , and  all  possible  configura- 
tions are  equally  probable.  The  viscous  forces  are  then  reduced  to 
zero.  So  far  as  the  individual  molecules  are  concerned,  this  state 
is  precisely  the  same  as  the  initial  state  even  though  the  macroscopic 
dimensions  of  the  material  have  changc-d  permanently. 

Let  us  consider  finally  what  happens  during  an  experiment  in 
stress  relaxation  in  which  the  elastomer  is  suddenly  extended  to  a 
constant  length.  If  the  extension  is  rapid  enough  relative  to  the  rate 
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of  diffusion,  the  bonds  are  distorted  and  new  degrees  of  freedom  are 
introduced  into  the  molecular  configurations.  In  this  state,  ^ cannot 
be  described  solely  in  terms  of  6,  X . After  the 

extension  is  fixed,  the  bonds  return  quickly  to  their  equilibrium  length, 
and  the  stress  simultaneously  decreases  from  a large  value  to  a more 
moderate  one.  At  this  moderate  stress,  the  molecular  configurations 
can  again  be  expressed  in  the  internal  coordinate  space  of 
Q,  *?  > X j...  X2f^(  , and  the  intermolecular  forces  are  entirely 

viscous  in  nature.  At  still  later  times,  the  molecules  slowly  diffuse 
back  to  the  equilibrium  sta*;e,  The  force  decays  to  zero  during  this 
process  in  the  same  way  as  it  decays  during  the  relaxation  described 
first  (constant  length  after  slow  elongation),  and  may  be  calculated 
from  the  solution  of  the  diffusion  equation. 

It  is  apparent  from  this  discussion  that  the  diffusion  theory 
of  Section  II  correctly  desc* *ues  the  stress  relaxation  only  at  fairly 
long  times  after  extension.  Unless  changes  in  the  bond  lengths  and  bond 
angles  are  introduced,  the  theory  cannot  be  used  to  describe  the  initial, 
rapid  decay  of  the  stress  [10], 


10,  A,  V,  Tobolsky,  J.  Am.  Chem.  Spc . 74,  3785  (1952)  points  out  that 
there  are  two  distinct  phenomena  occurring  during  the  stress 
relaxation. 
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IV.  GENERAL  SOLUTION  OF  THE  DIFFUSION  EQUATION 

In  order  to  find  a general  solution  for  s(t)  in  terms  of 
f(t)  , let  us  assume  that  e(t)  is  given  by  (the  real  part  of) 

i Out 

6 = €o  e • (36) 

The  (complex)  amplitude  £a  is  independent  of  the  time.  From  Equation  (14 ), 
• ^ 

the  potential  is  then 

V = (37) 


The  resistance  constants  will  all  be  equal  to  a constant  in 

a solution  so  dilute  that  there  are  no  interactions  between  the  molecul.  s. 
Whether  is  equal  to  or  not,  V may  be  written  as 

V-V.e1"** 


where 


v.-^i  |(pr-2»f). 


For  a potential  of  the  above  form,  f may  be  expanded  in 


Fourier  series  in  the  time: 

“ in«t 

4 = 2_  **  e » (40) 

r\  - Q 

where  fn  is  independent  of  the  time  but  is  a function  of 

XaH_,  • if  ea  is  small,  the  potential  of  any  possible 
configuration  will  be  small  compared  with  keT  . From  Equation  (19), 
it  follows  that  fh  is  much  smaller  than  fh_,  , and  the  higher  terms 

^The  subscript  on  the  potential  will  be  omitted  hereafter. 
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in  the  Fourier 

series  for  ^ may  be  neglected: 

(41) 

Substitution  of 

V and  f into  Equation  (19)  gives 

v-  D Vf0  = o 

(42) 

and 

V0-[Vf+foVl/AoT]  = 

(43) 

If  V 

were  identically  zero,  £ would  reduce  to  , 

Con- 

sequently,  f0  must  represent  the  state  of  equilibrium  and  be  independent 
of  ©,  i>  , Y-  2 j Xj  , X2b_(  . When  normalized, 

^°=Lfd£0  = \ / 2 ( 2rr)  , (44) 

and  is  identical  with  the  f0  introduce  previously. 

Let  us  now  define  the  orthogonal  eigenfunctions  of  the 
diffusion  operator  V D-V  : 

V-P-V^  = °‘  (45) 

The  )AX  are  a complete  set  of  all  the  functions  of 

satisfying  Equation  (45).  Each  eigenfunction  ^x  has  a corresponding 
eigenvalue  X . The  eigenfunctions  are  to  be  normalized  as  well  as 
orthogonal: 

(46) 

where  ^ is  the  complex  conjugate  of  and  ox  is  Kronecker's 


delta. 


28 


X>  may  be  expanded  in  terms  of  the  eigenfunctions: 


From  Equations 
independent  of 

by 


= 

(39)  and  (46) 

0,^*2, - 


2_  dx^*  • 


, the  expansion  coefficients  dx 
, Y-2n_,  as  well  as  the  time. 


(47) 

, which  are 
are  given 


(48) 


The  time  dependent  part 
expanded  in  terms  of  the 


of  the  configurational 
: 


X 


probability 


also  may  be 
(49) 


Substituting  the  expansions  of  V and  f,  into  Equation  (43),  using 
Equation  (45),  and  integrating  over  the  internal  coordinate  space  after 

• >i  * 

multiplication  by  , we  obtain 


< wf,  *?»  / fc  J 

4 k^T (l  + lw/X)  x 


(50) 


Finally,  substituting 

4 = f o 


bx  into  the  expression  for  $ 
«■  ^ c» g.  fu  y e^t- 

4 Ic.T  (|  -»-  iu./A) 


we  have 


By  Equations  (45)  and  (51), 

V-D-Vf  = 


V • D • Vjf  is 

i^.goiTc  y elwt 

4k0T  y-  (l+  ioj/X) 


(51) 


(52) 
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The  complete  conjugate  of  the  potential,  by  Equation  (47 ),  is 

* VA*uT  ^ 

V - ~zjt~  Ldv  n e • 


Introducing  these  last  two  expressions  into  Equation  (33)  together  with 
de/dt  from  Equation  (36)  and  using  the  orthonormality  of  the  , 


we  have 


/f\  — ^ nou>  *o  S-‘»  \ d v d* 

•5(i)  = - (,  + ;w7r) 


This  equation  gives  the  amplitude  and  phase  of  the  stress  that  must  be 
applied  to  a cylinder  of  the  elastomer  in  order  to  produce  a sinusoidal 
extension  of  amplitude  . 

To  the  approximation  that  the  stress  is  a iinear  function  of 
the  strain,  a convenient  way  to  express  the  general  experimental  relation- 
ship between  the  stress  and  the  extension  is  in  terms  of  a distribution 
of  relaxation  times  E'M  , c'(t)  is  defined  by 


f*  ( -(t  - t')/r  , d i(t')  , 

set)  = | e L(T)-r-dtdr. 

v/  J 


= -oo 


. tU/l 

"When  e(t)  is  given  by  the  real  part  of  cQ e , the  succeeding 

equations  for  V,  f , and  s(t)  should  also  be  prefaced  on  their  right 
sides  by  the  words  "the  real  part  of."  Equation  (33)  for  the  stress 
should  be  written  as  the  real  part  of  V*  , multiplied  by  the  real 
part  oi  5/.  D-V^f  , and  divided  by  the  real  part  of  de/dt  . How- 
ever, V is' always  in  phase  with  de/dt  , and  "the  real  part  of" 
may  be  omitted  in  all  the  equations.  It  should  also  be  noted  that 
Vfc  is  determined  by  changing  the  factor  in  V dependent  upon  the 
configuration  of  the  molecule  and  does  not  involve  a change  in  the 
complex  dependence  upon  t . 


l 
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The  sinusoidal  solution  of  the  diffusion  equation  given  by  Equation  (54) 
is  sufficient  to  calculate  the  distribution  of  relaxation  times.  Let  us 
define  /2(X)  as  the  sum  of  all  the  terms  d,-  having  X'  less 

than  X : 


X't  X 


(56) 


In  practice  it  is  nearly  impossible  to  perform  this  summation  exactly, 
and  O(X)  must  be  approximated  by  a continuous  function  of  X . To 
this  approximation,  its  derivative  exists: 


H(M  = 


dO  (x) 
d X 


(57) 


By  introducing  H(M  into  Equation  (54),  we  may  express  s(t)  as  an 
integral  instead  of  a summation: 


s(t) 

Substituting  a new  variable 
into  this  equation,  we  have 


— —— — — I — e 


lb*.  I j (l  +•  «.t a/X) 

W.  5 


r s T 


(53) 


(59) 


,Tt  = inB<0C.g.a  ^ 

s It  k T J r3(ir  loir) 


For  sinusoidal  motion.  Equation  (55)  becomes 


(60) 


s(t)  - 


- tw  e 


■ f T E'fr)  dr 

° J ( I + ;wt)  e • 

rzo 


(61) 


I 
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If  Equations  (60)  and  (61)  are  to  give  the  same  amplitude  and  phase  for 
the  stress,  EL  (t)  must  be  identical  with 


L\r)  = 


n.  H(?) 
l!»k#T  r4  ' 


(62) 


This  equation  completes  the  derivation  of  the  general  theoret- 
ical relationship  between  the  stress  and  the  strain,.  The  method  for 
finding  E(t)  from  the  diffusion  tensor  D may  be  summarized  as 
follows:  (1)  the  eigenfunctions  of  V*Q-*7  are  calculated  from 

Equation  (45);  (2)  from  the  , the  expansion  coefficients  of 

the  potential  are  obtained  from  Equation  (48);  (3)  the  function  n«./W 
is  calculated  from  the  dx  by  Equations  (56)  and  (57);  and  (4)  H(M 
is  substituted  into  Equation  (62)  to  find  f(r)  . The  solution  for 
ELC"*")  may  then  be  substituted  into  Equation  (55)  to  find  the  stress 
corresponding  to  any  form  of  e(tr)  , whether  sinusoidal  or  not. 


V.  SPECIFIC  SOLUTION  FOR  THE  STRESS  RELAXATION 


In  this  section,  the  distribution  of  relaxation  times  is 
computed  from  the  diffusion  tensor  obtained  originally  by  Kirkwood  and 
Fuoss,  The  stress  relaxation  at  constant  extension  is  then  calculated 
from  £("*”)  and  compared  with  published  experimental  results  for  poly- 
isobutylene . 

The  first  step  in  finding  £(Y)  from  the  general  derivation 
in  Section  IV  is  to  evaluate  D,  the  diffusion  tensor.  D is  defined 
as  k0T  p ' , where  /O  expresses  the  relationship  between  the  viscous 

torques  and  the  rates  of  change  of  the  angular  coordinates 
6 ,4*  ^2n-l  : 

Ti='Le^i-  (63> 

j 

The  vesistance  tensor  p has  been  evaluated  by  Kirkwood  and  Fuoss  ; 
their  derivation  is  reproduced  in  Appendix  I.  The  tensor  p is  dependen 
upon  the  resistance  constants  of  the  individual  elements  of  the 

chain.  It  is  assumed  in  Appendix  I that  all  of  the  ar8  equal  to 

, where  (introduced  and  defined  approximately  on  page  25)  is 

the  resistance  constant  of  an  element  moving  in  a liquid  composed  of 
the  unpolymerized  elements.  In  other  words,  it  is  assumed  that  the 
viscous  force  on  a polymer  chain  is  the  same  whether  or  not  the  sur- 
rounding fluid  is  polymerized. 

The  resistance  tensor  p resulting  from  this  assumption  is 
a complicated,  nondiagonal  function  of  the  angles  defining  the 


*See  Appendix  I of  [7]. 
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configuration  of  the  chain.  It  is  so  complicated  that  exact  evaluation 
of  its  inverse  is  impracticable.  However,  we  cam  approximate  £ by 
noting  that  is  much  smaller  than  in  Equation  (41).  In  other 

words,  the  distribution  of  molecular  configurations  is  very  close  to 
the  equilibrium  distribution  . Thus,  there  should  be  only  a small 
change  in  E(r)  if  yO  is  replaced  by  , its  average  over  all  pos- 
sible configurations  of  the  molecule.  To  this  approximation,  the 
diffusion  tensor  becomes 

D 2 kT (£)’.'  (64) 

Evaluation  of  D is  quite  simple  because  /5  is  diagonal  in  the  internal 
coordinate  space.  From  Equations  (130)  and  (131)  in  Appendix  I,  the 
nonzero  components  of  0 are 

p8a=5^  = D0/2  (65) 

and  D,.  = D.[  I - 1 ~|  ] , (66  ) 

where  Q,  “ (67) 


a,  is  the  bond  length  and  n is  half  the  number  of  elements  per  mole- 
cule. (For  substituted  polyethylenes,  n is  equal  to  the  number  of 
monomers  per  molecule.) 

From  Equations  (133)  and  (134),  the  corresponding  eigen- 
functions *k  of  the  diffusion  operator  are 

IU  r(2*-H)(i-Kl)!  1^/0  )n  + l pm»/  i"\Xt 

^ = i 4-n-(iTi  «.orJ  k2irj  ^(cosd)<?  Ue 


(68) 
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and  the  eigenvalues  are 

X=  Da{riU+0  + ]T  [l-|1^L|jmi  }•  (69) 

The  are  the  associated  Legendre  polynomials.  The  toteil  set  of 

eigenfunctions  includes  the  ^ for  all  positive  A,  all  ma  whose 
absolute  values  are  less  than  or  equal  to  & , and  all  positive  and 
negative  tni  , The  above  results  are  those  obtained  by  Kirkwood  and 
Fuoss  [7], 

The  next  problem  is  to  find  the  expansion  coefficients  dx 
of  the  potential.  From  Equation  (48),  the  c!^  are  equal  to 

dx  = ■ Z ~ (70) 

when  . We  might  note  that  the  potential  is  a quadratic 

function  of  the  position  of  the  elements.  This  case  contrasts  with 
that  for  an  electric  field,  where  V is  linearly  dependent  upon  the 
positions  of  the  polar  side  groups.  Because  the  potential  is  quadratic, 
evaluation  of  the  dx  is  quite  laborious.  A summary  of  the  compu- 
tations for  the  dx  is  given  in  Appendix  II. 

After  the  have  been  determined,  the  next  step  is  to 

find 

T (7!) 

A (A)  is  the  sum  of  the  squares  of  all  dx.  whose  corresponding 
eigenvalues  are  less  than  A . It  is  nearly  impossible  to  evaluate 
this  summation  exactly,  but  it  can  be  approximated  in  the  following 
way.  Let  us  consider  all  the  dx  having  some  particular  absolute 
value.  These  dx  have  corresponding  eigenvalues  (given  by  Equation  (69)) 
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that  are  all  nearly  equal.  It  la  not  too  much  in  error,  then,  to  assume 
that  the  dA  having  the  same  absolute  value  also  have  the  same  eigen- 
value. To  this  approximation,  is  given  by 


UW  = [ N(SHd*i  . <72> 

where  N(M  is  the  number  of  dx  having  some  particular  absolute 
value.  The  corresponding  average  eigenvalue  is  X . This  approximation 
to  ncx)  is  similar  to  one  used  by  Kirkwood  itnd  Fuoss  in  their  treat- 
ment of  the  dielectric  dispersion. 

is  evaluated  in  Appendix  III.  From  Equations  (187) 


and  (189),  we  have 


•)  (c*->r 


}• 


(73) 


where  c = 3A/2.fiD0.  (74) 

Onco  is  known,  ri(A)  and  E(r)  can  be  found  from 

Equations  (57)  and  (62).  The  distribution  of  relaxation  times  is* 


r'M  - 8N0(lc0T)V 

J~^.Q?(c  + 3)s 


(75) 


When  expressed  in  terms  of  t , the  variable  c is 

C = T./'r,  (76) 

where  % = C.an/2  k.T'  (77) 

No  is  half  the  total  number  of  eleiments  in  a unit  volume  of  the  polymer, 
and  is  independent  f the  molecular  v?,i(;ht.  (For  substituted  polyethylenes. 


^Equations  (75),  (78),  and  (79)  are  correct  only  if  i;/n  « t 
If  T < T,/n  or  r > T.f-ii  E'(r)  is  zero.  Because  ri  is  very  large, 
these  limits  can  be  neglected  in  the  calculations  for  s(t). 
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N0  is  the  number  of  monomers  per  unit  volume „)  The  time  Ta  is  a 
function  of  both  the  temperature  and  the  molecular  weight. 

The  distribution  of  relaxation  times  is  plotted  versus  r/Y0 
in  Figure  2,  with  logarithmic  scales  for  both  variables.  The  limiting 
value  of  E(t)  for  small  r is 

£ 

E'(T)=W  . (78) 

3 s0a 


As  r increases  from  zero,  E'(t)  rises  slowly  to  a maximum  at  TaO.lT,, 
and  then  drops  off  rapidly.  When  T is  large 


(79) 


The  next  problem  is  to  determine  how  closely  the  theoretical 
distribution  describes  the  exper iment al  properties  of  noncross-iinked 
elastomers.  We  shall  calculate  the  stress  relaxation  from  E.'(t)  «nd 
compare  .ut  with  the  relaxation  of  a typical  elastomer.  One  of  the  most 
nearly  complete  experimental  investigations  [11]  has  been  of  polyiso- 
butylene, and  wo  shall  compare  the  theoretical  results  with  the  data 
for  this  polymer. 

From  Equation  (55),  the  stress  relaxation  s(i)  at  constant 
extension  € is  given  by 

co 

s(t)|  = ej  e E'(T)dr. 


T=o 


(80) 


11,  R.  S.  Marvin,  "Interim  Report  on  the  Cooperative  Program  on  Dynamic 
Testing,"  National  Bureau  of  Standards,  1951;  and  E,  R.  Fitzgerald, 
L,  D,  Grandine,  Jr.,  and  J.  D.  Ferry,  J,  Appl.  Phys . 24 . 650  (1953), 
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The  time  t .is  measured  from  the  ina.anv  j»f  extension.  To  find  5 6^ 
it  is  necessary  to  substitute  E(t)  into  this  equation  and  integrate 
over  r for  various  values  of  t . 

The  calculations  for  s(.t)  are  described  in  Appendix  IV,  ana 
s(t)  is  plotted  in  Figure  3.  In  the  figure,  the  stress  is  normalized 
to  unit  force  at  t-O  by  dividing  by 

s(°)|€=  yN.nkJc.  (81) 

No  explicit  equation  for  can  be  given  because  the  integration 

is  performed  graphically. 

Also  shown  in  Figure  3 is  the  composite  experimental  stress 

relaxation  of  an  unfractionated  polyisobutylene  at  30°C  [12],  The 

average  molecular  weight  is  6,600,000.  The  curve  has  been  normalized 

6 2 

by  dividing  the  original  data  by  10.4  x 10  dynes/cm  -unit  extension. 

The  time  scale  for  the  theoretical  stress  relaxation  has  been  chosen 
so  that  the  two  curves  coincide  at  a normalized  stress  of  0.50 
( T0  = 102’04  hours). 

The  experimental  stress  relaxation  may  be  divided  into  two 

quite  distinct  regions.  One  region  includes  the  slow  relaxation  of 

7 2 

stresses  less  than  10  dynes /cm  -unit  extension.  This  relaxation  is 
due  to  the  molecular  diffusi:  •>  treated  ir.  this  thesis.  The  other  region 
includes  the  stresses  larger  than  10  dynes/cm  -unit  extension.  These 
forces  decay  very  rapidly  at  30°C.  and  have  almost  entirely  disappeared 


12,  The  composite  curve  has  been  formed  by  superposition  of  the  data 
plotted  in  Figure  2 of  R.  D.  Andrews  and  A.  V.  Tobolsky, 

£.  Polymer  Sci.  7,  221  (1S51). 


TIME  IN  HOUFtS 
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10  hours  after  extension.  For  polyisobutylene,  this  relaxation  is 
not  experimentally  observable  unless  tho  tempera+ure  is  much  lower  than 
30°C,  No  attempt  has  been  made  to  include  the  more  rapid  relaxation 
in  the  theoretical  discussion,  and  th-  calculated  curve  may  be  compared 
only  with  the  relaxation  at  relatively  large  times.  For  this  reason, 
the  normalization  factor  has  been  chosen  so  that  the  stress  due  to 
diffusion  is  unity  when  t is  equal  to  zero. 

It  is  apparent  from  the  figure  that  the  stress  decays  much 
more  slowly  than  is  predicted  theoretically.  This  result  implies  that 
the  relaxation  times  given  by  Equation  (75)  are  not  distributed  in  T 
as  widely  as  would  be  necessary  to  describe  the  experimental  data  for 
polyisobutylene. 

The  theoretical  and  experimental  stress  relaxations  also  do 

not  agree  in  their  dependence  upon  the  molecular  weight.  Experimentally, 

3.3 

the  time  scale  of  the  relaxation  is  proportional  to  H [12,13],  The 
magnitude  of  the  force  is  not  dependent  upon  PI.  In  other  words,  a 
plot  of  s(t)  versus  "t/  M is  independent  of  the  molecular  weight. 
From  Equations  (77)  and  (81),  it  may  be  seen  that  the  calculated  results 
are  quite  different:  both  the  force  and  the  time  scale  are  proportional 

to  the  first  power  of  M.  A plot  of  the  calculated  s(t)/M  versus  "fc/PI 
is  independent  o:  the  molecular  weight. 

Ko  are  more  successful  in  comparing  the  predicted  and  exper- 
imental dependence  cn  temperature.  The  only  parameter  affecting  E(t) 


13,  T,,  G.  Fox,  Jr.,  anu  P.  J.  Kiory,  J_,  Am,  Chora,  See.  70,  2384  (1948), 
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that  can  depend  upcn  T (other  them  the  temperature  itself)  is  the 

resistance  constant  • Whether  or  not  the  elements  along  the  polymer 

chain  are  spheres,  is  proportioned  to  l<,  , the  viscosity  of  the 

A 

liquid  composed  of  unpolymer lzed  elements  . For  nearly  all  unpolymer - 
ired  liquids,  the  viscosity  is  exponentially  dependent  upon  l/T,  and 


presumably  varies  in  the  same  way  here: 

. A/ftT 

£ * * e 


(82) 


If  this  temperature  dependence  is  substituted  into  Equation  (77),  we 
see  that 


r * 


-l  A/KT 
’ e 


(83) 


Also,  from  Equation  (81)  for  s(°)  , the  calculated  stress  is  directly 

proportional  to  the  absolute  temperature  and  independent  of  C>  . The 

ft  j rt 

proportionality  of  s(o)  to  T ar.d  of  X to  a ' correspond  to 
the  temperature  dependence  observed  experimentally  for  polyisobutylene 
[11,12], 

It  is  possible  to  find  the  dependence  of  the  stress  upon  the 
temperature  without  calculating  s (t)  explicitly,  -fl(\)  consists 
of  a partial  sum  of  the  expansion  coefficients  of  the  potential. 

Change  of  the  temperature  will  alter  the  relaxation  time  i/a  of  each 
expansion  coefficient  but  will  not  change  its  magnitude.  This  effect 
implies  that  J"1  is  equal  to  some  function  of  Xx  and  not  just  A 
alone ; 

n =n(Xx).  (84) 


Tor  a discussion  of  the  relationship  between  the  viscosity  and  the 
resistance  constants,  see  [6,8], 
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I 

I 

{ 


From  Equations  (57),  (62),  (80),  and  (84),  ve  obtain 


s(t) 


f = 


where  the  variable  £ denotes  t/t#  , and 


Q(P)  = p 


r4dA(W 


dfi 


(85) 


(86) 


From  Equations  (83)  and  (85),  it  follows  that  a plot  of  s(-fc)/T  versus 

, -T-  -AiR T _ . 

t I 4 is  independent  of  the  temperature.  The  stress  at  t —O 

is  therefore  proportional  tu  T , and  the  time  scale  of  the  relaxation 
depends  exponentially  upon  l/T  . This  temperature  dependence  satisfies 
the  usual  reduction  method  [14]  for  relating  the  mechanical  properties 
of  elastomers  at  different  temperatures,  and  agrees  with  the  experimental 
results  for  polyisobutylene. 

From  the  above  discussion,  it  may  be  seen  that  the  dependence 
of  s(t)  upon  ths  temperature  is  not  related  to  the  specific  form  of 
the  diffusion  tensor.  On  the  other  hand,  the  dependence  of  s('t)  upon 
the  time  and  the  molecular  weight  does  vary  with  jg  . It  would  there- 
fore appear  that  the  general  solution  for  £(1)  in  Section  IV  may  be 
correct,  but  that  the  approximations  in  this  section  to  find  g are  in 
error. 

An  important  source  of  error  is  the  assumption  that  the 
resistance  tensor  of  a molecule  is  the  same  whether  or  not  the  surrounding 


14,  R.  S,  Marvin,  E,  R.  Fitzgerald,  and  J.  D.  F'erry,  _J.  Appl.  Phys.  21. 
197  (1950);  and  J.  D.  Ferry,  E.  R,  Fitzgerald,  M,  F,  Johnson,  and 
L,  0,  Grandine,  Jr,,  J.  Appl.  Phys . 22.  717  (1951);  see  also 
F,  Scnwarzl,  and  A.  J,  Staverraan, _J.  Appl.  Phys.  23,  838  (1952). 


fluid  is  polymerized.  Furthermore,  the  approximation  of  O to  0 may 
be  inaccurate  enough  to  change  £ (t)  appreciably.  Both  of  these 
approximations  tend  to  make  £(T)  too  narrow  a function  of  T , From 
this  point  of  view,  the  theoretical  stress  relaxation  in  Figure  3 may 
be  considered  as  an  upper  limit  upon  the  rate  at  which  the  relaxation 
can  proceed. 

It  might  be  pointed  out  that  the  theoretical  distribution  of 
relaxation  times  should  agree  much  better  with  the  extension  properties 
of  a dilute  solution  of  a polymer.  In  a solution,  the  fluid  surrounding 
each  polymer  molecule  is  unpolymerized  and  is  equal  to  . 
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V I . CORRECTION  OF  THE  MECHANICAL  PROPERTIES  FOR  CHAIN  ENT ANGLEMENTS 


From  the  previous  section  it  may  be  concluded  that  the  cal- 
culated stress  relaxation  does  not  agree  with  the  experimental  results 
for  a typical  elastomer,  and  that  the  primary  reasons  for  this  discrepancy 
are  the  approximations  used  to  find  the  diffusion  tensor.  To  obtain  a 
more  precise  f('r')  , it  is  necessary  to  evaluate  the  effects  of 
polymerization  of  the  surrounding  molecules  upon  the  resistance  tensor 
p of  a polymer  chain. 

Let  us  examine  qualitatively  the  interactions  of  the  chains 

by  considering  a rotation  around  some  azimuth  angle  X,  , If  all  the 

J 

other  angles  along  the  same  chain  are  fixed,  one  part  of  the  molecule 
(on  that  side  of  the  j-th  bond  containing  the  center  of  the  chain) 
remains  stationary.  The  other  and  smaller  half  rotates  around  the  bond 
as  a rigid,  irregular  rod.  If  the  surrounding  molecules  are  small  (as 
in  a dilute  polymer  solution),  the  forces  on  the  chain  can  bo  calculated 
from  the  resistance  constant  f?m  of  each  element.  Evaluation  of  these 
forces  leads  directly  to  the  diffusion  tensor  used  in  the  previous 
section. 


In  the  unp^ •A.vt.lcizcd  polymer,  however,  the  molecules  are 
entangled  with  each  other.  The  part  of  the  chain  rotating  around  the 
j-th  bond  carries  with  it  a large  number  of  surrounding  elements,  and 
the  resistance  to  the  rotation  is  much  larger  than  the  corresponding 
force  in  the  unpolymer ized  fluid.  If  the  effective  resistance  constant 
of  each  element  includes  the  viscous  forces  on  all  the  chains  wound 
, is  much  larger  than  £ . Furthermore,  if 


around  that  element 
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n , the  number  of  chains  caught  on  the  i-th  element  is  larger  than 
the  number  caught  on  the  (i-l)th  element.  This  comes  about  because  the 
chains  are  free  to  slip  off  the  end  of  the  rotating  molecule.  As  a 
result,  the  effective  C;  are  larger  near  the  center  of  the  chain  than 
they  are  at  the  ends. 

He  might  summarize  this  discussion  by  suggesting  that  the 
diffusion  tensor  may  be  calculated  in  the  same  way  as  in  Appendix  I, 
but  with  4*.  introduced  as  a function  of  i . 

We  shall  not  attempt  to  determine  the  effect  of  the  chain 
entanglements  quantitatively.  Instead,  we  shall  assume  some  ad  hoc 
dependence  cf  upon  i-  and  calculate  the  corresponding  stress  relax- 

ation. The  purpose  of  this  procedure  is  to  find  a simple  set  of 
resistance  constants  correctly  describing  the  mechanical  properties  of 
elastomers.  Once  the  proper  are  determined,  they  should  be  appli- 

cable to  the  solution  of  the  diffusion  equation  for  all  types  of  external 
forces,  whether  mechanical  or  not. 

Let  us  assume  that  is  given  by 

Kalp 

^ KA12-1-/,  (87) 

where  A and  p are  constants  independent  of  the  temperature.  These 
diffusion  constants  are  symmetric  around  the  center  of  the  chain, 
increase  rapidly  toward  its  center,  and  are  a multiple  A of  ^ at 

O 

the  ends  of  the  molecule.  From  Equation  (14)  it  is  apparent  that  the 
potential  as  well  as  the  diffusion  tensor  is  modified  by  this  assump- 


tion. 


J 
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The  new  D , dx , and  fi(X)  are  calculated  in  Appendix  V. 
If  the  function  £(c)  is  defined  by 

..  ( ) 'r*C  (LtL)  litll) 

S(c)=  (C43)1  j 


nex)  may  be  expressed  as 


iP*4-  4 A2 

linr  a A 


n(x)  = 

where  c is  equal  to  V / 7 and 

r j =»  C.Aan5'  /(p^2)HJ. 

The  corresponding  distribution  of  relaxation  times  is1 

dH(c) 


zb)  - 


N.(p-^2)3(I<,t/c4 


3(p*i)3(2p4  3)'?.A£^  <JC  * 


(88) 


(89) 


(90) 


(91) 


If  T is  very  small  compared  with  t.'  , 


Ma(?42)3(k0Tf  ( T \{^]  , 

ElT)-3(p+.)3«>  r<<° 


(92) 


If  T is  much  larger  than  t0'  , 


rV  ei  A/0 (p-^2)  (k,T)  ( r \~J 

<?(p-Kt)4^Aa  np  ' t; /• 


r >? 


(93) 


liquations  (91),  (92),  and  (93)  are  correct  only  if  \/hf’*  « ~r  «^n- 
When  t < 'rl/n^''  or  t ■»  T.'n  ^ E.'(T)  is  zero.  Tteie  limits 

can  be  neglected  in  the  calculations  for  s(t:)  . 


i 
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From  Equation  (80) , the  stress  relaxation  is 


S 


(p*2)*M  = ni<TC  r -(- ^.)c  2 dfi(c)  jc 
3(rl)3(2p+3)  ) * dc 


(94) 


Two  important  results  may  be  obtained  directly  from  these 
equations.  In  the  first  place,  the  temperature  dependence  of  the  stress 
relaxation  is  the  same  as  in  the  previous  section.  The  stress  at  small 
times  .is  proportional  to  the  absolute  temperature,  and  the  time  scale  is 
dependent  exponentially  upon  \/t  . Secondly,  V is  dependent  upon 
the  molecular  weight  raised  to  the  prl  power,  provided  that  A and  p 
are  independent  of  M.  The  form  of  this  dependence  corresponds  quite 
well  to  the  experimental  relationship  that  the  time  scale  of  the 
relaxation  is  proportional  to  fi  ' . 

The  stress  relaxation  of  polyisobutylene  may  be  fitted  quite 
accurately  by  setting  p equal  to  five.  The  calculations  for  ECv)  and 
s(t)  when  p is  five  are  described  in  Appendix  VI,  From  Equation  (220), 


Ur)  = 


0.04072  N.(KT)2c4 [ 

e.aJnJAC-.3)'0'n 


7-287 +-(c+l)fe[^(6 


3 10  Q-3b\  /c+3  XT 1 ) 

(c+l)  (c  + l)1/  \c  -W  / JJ* 


(95) 


This  distribution  of  relaxation  times  is  plotted  against  the  logarithm 
of  t / vo  in  Figure  4.  The  corresponding  stress  relaxation  is  plotted 
in  Figure  5,  with  the  force  at  t = 0 normalized  to  unity  by  dividing  by 


s(0)|  =0.3  4 8 N#nl«.Te. 


(96) 


Also  shown  in  this  figure  is  the  experimental  stress  relaxation  plotted 
in  Figure  4 (polyisobutylene  for  M = 6,600,000  and  T=  30cC)  [12], 

The  experimental  curve  is  again  normalized  by  dividing  by 


TIME  IN  HOURS 


46 


6 2 

10,4  x 10  dynes/cm  -unit  extension;  ~rj  is  chosen  so  that  ine  curves 
coincide  at  a normalized  stress  of  0.50  ( \ = 104*^  hours). 

The  theoretical  curve  agrees  with  the  composite  exper iment al 
curve  to  within  0.05  over  nine  decades  of  time.  Only  at  very  small 
times  after  extension  do  the  two  curves  disagree.  At  these  times,  the 
relaxation  is  not  produced  by  diffusion  of  the  azimuth  angles  of  the 
molecules.  It  may  be  concluded  that  a set  of  diffusion  constants  givon 
by  the  rule 


KAi5 


i.  * o 
L ? f * I 


(97) 


predicts  a stress  relaxation  in  excellent  agreement  with  the  properties 
of  polyisobutylene.  The  temperature  dependence  is  also  correctly 
described  py  this  law. 

the  molecular  waight  dependence  predicted  by  these  does 

not  quite  agree  with  the  experimental  results.  In  the  first  place,  the 
calculated  V is  proportional  to  the  sixth  power  of  M,  rather  than 
the  3,3  power  observed  experimentally.  This  disagreement  is  not  totally 
unexpected,  however.  We  have  selected  a value  of  p which  corrects  the 
time  dependence  not  only  for  the  effect  of  the  chain  entanglements  but 
also  for  the  substitution  of  D for  0 , The  latter  approximation 
changes  the  time  dependence  of  s ( t)  without  altering  the  dependence 

upon  the  molecular  weight.  Consequently,  a p of  five  overcorrects 
for  M . If  A were  independent  of  H , we  might  conclude  that  about 


*The  introduction  of  p also  corrects  for  the  approximation  to  f)(X) 
in  Equation  (72),  and  for  the  distribution  of  molecular  weights  in 
the  unfractionated  polymer.  Both  of  these  factors  change  the  shape 
of  s(t)  without  altering  the  effect  of  M , and  therefor  s tend  tc 
make  the  predicted  dependence  of  >0'  upon  M too  large. 
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half  of  p corrects  for  the  effects  of  chain  entanglements , while  the 
other  half  corrects  for  the  approximation  to  D . A more  accurate 
solution  of  the  diffusion  equation  than  is  given  here  should  lower  the 
dependence  of  upon  M , and  agree  more  closely  with  the  experimental 

results. 

The  calculated  stress  relaxation  depends  in  a second  way  upon 
the  molecular  weight.  The  magnitude  of  s(o)  from  Equation  (96)  is 
proportional  to  (1  , whereas  experimentally  it  does  not  vary  with  fi  , 

One  might  suspect  that  s(o)  is  independent  of  M for  an  unplasticized 
polymer  because  the  chain  entanglements  act  as  temporary  cross-links 
between  the  molecules.  In  other  words,  the  rapid  motion  of  one  part 
of  a chain  may  be  entirely  divorced  from  the  rapid  motion  of  another 
part  of  the  same  chain.  Tho  stress  at  small  times  would  then  be  inde- 
pendent of  the  total  length  of  the  molecule.  Only  over  relatively  long 
periods  of  time  would  the  motions  of  the  various  parts  of  a chain  be 
related.  During  such  motion,  the  sole  effect  of  the  entanglements 
would  be  to  change  the  resistance  constants  of  the  elements.  Since 
the  stress  relaxation  is  primarily  a slow  phenomenon,  the  treatment 
given  here  may  be  essentially  correct  with  regard  to  the  shape  of  the 
relaxation,  even  tr.ough  it  is  not  correct  as  a calculation  for  s(o)  . 

It  might  be  pointed  out  again  that  the  use  of  E^r)  is  not 
restricted  +o  a calculation  of  the  stress  relaxation.  Equation  (55) 
gives  s(t)  corresponding  to  almost  any  e(t)  whatsc- ^er.  For  example, 
we  may  use  E,('>')  to  compute  the  stress  required  for  a constant  rate 
of  extension.  The  only  limitation  upon  Equation  (55)  is  that  de/dt 
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must  not  bo  too  large.  In  goner al,  the  published  data  for  polyisobutylene 
are  internally  consistent,  and  the  E (~)  which  gives  the  proper  stress 
relaxation  also  will  predict  the  correct  response  to  other  types  of 


mot  ion 


49 


VII.  CORRECTION  OF  THE  DIELECTRIC  DISPERSION  FOR  CHAIN  ENTANGLEMENTS 

The  diffusion  equation  for  £ (Equation  20)  can  give  the 
response  of  an  elastomer  to  forces  other  than  those  accompanying 
mechanical  extension.  The  only  restriction  upon  the  equation  is  that 
the  molecular  motion  must  consist  of  diffusion  in  the  angles 
0 , ^ , X,  , X . The  diffusion  equation  is  applicable  to  the 

dielectric  dispersion  of  polar  polymers,  for  example,  if  the  potential 
is  changed  to 

k Lklt 

V = £ • « • (98) 

Fe  is  the  amplitude  of  the  local  electric  field  of  frequency  u>/2tt 
and  p is  the  total  dipole  moment  of  the  polymer  molecule. 

The  diffusion  equation  for  this  potential  has  been  solved 
by  Kirkwood  and  Fuoss  for  a polymer  in  which  the  dipoles  are  attached 
rigidly  to  the  atoms  making  up  the  molecular  chain.  Polyvinyl  chloride 
is  an  example  of  such  a polymer.  Their  calculation  is  intended  to 
apply  to  a solution  dilute  enough  so  that  all  chain-to-chain  entangle- 
ments may  be  neglected.  Consequently,  they  use  the  diffusion  tensor 
derived  upon  the  assumption  that  all  the  resistance  constants  are  equal. 

The  substitution  of  D for  D and  the  approximation  to  -Q(\)  given 
by  Equation  (72)  are  also  employed  in  their  calculation. 

The  results  of  Kirkwood  and  Fuoss  may  be  expressed  as  a distri- 
bution of  electrical  relaxation  times  Q(t)  , For  a perfectly  fractionated 
polymer,  this  distribution  is  given  by 

-r» 


(99) 
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The  reduced  dielectric  dispersion  H(u>)  at  frequency  ^flir  is  the 
imaginary  part  of 


-OW= 

J t IWTI 


'r-o 


(100) 


3y  substituting  Q(t)  into  this  integral,  they  obtain 


- x*+irx-l],  (101) 

V * * / 

The  variable  x is  equal  to  , where  rra  is  the  relaxation  time 

defined  by  Equation  (77). 

H(<*>)  , divided  by  its  maximum  value  at  X = 1,  is  plotted 
against  the  frequency  in  Figure  6.  Also  plotted  in  this  figure  is  the 
audio  frequency  dispersion  of  unplasticized  polyvinyl  chloride  [15], 

The  data  were  taken  at  abcut  1C0°C.  The  frequency  scale  of  the  theoret- 
ical curve  has  been  chosen  so  that  Hmo„  coincides  with  the  maximum 

-2  99 

of  the  experimental  curve  ( % = 10  ’ seconds). 

The  experimental  dispersion  exhibits  a very  broad  maximum 
not  fitted  by  the  theoretical  curve.  We  may  conclude  here,  precisely 
as  for  the  stress  relaxation,  that  the  calculated  distribution  of 
relaxation  times  is  not  broad  enough  to  correspond  to  the  experimental 
data. 

It  is  not  surprising  that  the  results  of  Kirkwood  and  Fuoss 
disagree  with  the  experimental  dispersion  of  an  unplasticized  polymer, 
because  their  theory  applies  only  to  a dilute  solution  of  a polymer. 


15,  The  composite  experimental  curve  has  been  formed  by  superposition 
of  the  data  plotted  in  Figure  4 of  R.  M.  Fuoss, _J.  Am,  Chem.  Soc. 
63,  369  (1941), 


KIRKWOOD 
0 FUOSS  (p=0) 


51 


From  Section  VI,  however,  we  have  available  a set  of  resistance  constants, 
A f ) that  describes  the  effect  of  the  molecular  entanglements 
during  mechanical  extension  of  an  unplasticized  polymer.  Since  the 
same  entanglements  occur  during  the  application  of  an  electric  field, 
introduction  of  these  resistance  constants  into  the  Kirkwood  and  Fuoss 
theory  should  give  the  correct  dielectric  dispersion. 

The  solution  of  the  diffusion  equation  for  an  electric  field 
has  been  recalculated  using  the  modified  resistance  constants.  A summary 
of  the  computations  is  given  in  Appendix  VII.  The  new  distribution  of 
electrical  relaxation  times  is 

° ■ <102) 

where  is  given  by  Equation  (90).  H (co)  , divided  by  HmaK  , is 

plotted  in  Figure  6 against  the  frequency,  with  set  equnl  to 

10-1*99  seconds  so  that  Hma)L  coincides  with  the  maximum  of  the  exper- 
imental curve.  The  corrected  dielectric  dispersion  is  much  closer  to 
the  experimental  result  than  is  the  uncorrected  curve. 

It  might  be  pointed  out  that  the  dispersion  for  p equal  to 
five  is  not  symmetric  about  its  maximum.  This  suggests  that  the  Kirkwood 
and  Fuoss  theory  has  not  been  corrected  properly  for  frequencies  less 
than  2tt/t0'  . The  stress  relaxation  at  corresponding  times  is  too 

small  to  measure  accurately,  and  the  correction  is  not  applicable  in 

^Equation  (99)  is  correct  only  if  %/r'  <*  t <<  . Equation  (102) 

is  correct  only  if  %'/  h P* ' «.  "*•  <<  rj  n , These  limits  may  be 
neglected  in  the  calculations  for  H(u>). 
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this  region.  Over  the  three  decades  of  frequency  larger  them 
the  theoretical  curve  is  never-  more  than  0,08  below  the  experimental 
H(U))  • This  precision  is  as  good  as  can  be  expected  because  the  value 
of  p w^s  determined  from  the  mechanical  properties  of  a different 
polymer,  A more  accurate  check  or.  the  theory  would  be  the  calculation 
of  the  dielectric  dispersion  from  the  stress  relaxation  of  polyvinyl 
chloride  itself.  With  this  reservation,  the  modified  theory  is  an 
accurate  representation  of  the  electrical  properties  of  polyvinyl 
chloride, 

One  way  of  summarizing  the  dependence  of  H(uj)  upon  frequency 
is  to  gave  its  maximum  value,  because  decreases  as  the  width  of 

the  dispersion  increases.  For  polyvinyl  chloride,  is  about 

0,10  at  iOO°C.  In  the  Kirkwood  and  Fuoss  theory,  is  0,285  (for 

a fractionated  polymer).  When  the  correction  for  the  chain  entangle- 
ments is  -introduced  into  the  theoretical  derivation,  is  reduced 

to  0.143,  a value  almost  as  low  as  the  experimental  figure. 

The  predicted  temperature  dependence  of  H(oj)  also  agrees 
with  the  experimental  results.  For  both  of  the  theories  and  far 
polyvinyl  chloride,  the  frequency  of  maximum  dispersion  is  dependent 
exponentially  upon  \j T . T'ne  shape  and  magnitude  of  H {**>)  are 
independont  of  T.  Actually,  for  polyvinyl  chloride  there  is  a slight 
change  is  as  the  temperature  i=  raised,  but  this  effect  may 

be  due  to  other  types  of  dispersion  occurring  primarily  at  higher 
frequencies , 

From  Equation  (90)  we  can  show  readily  tha  calculated 

, .p+i 

frequency  of  the  maximum  dispersion  is  proportional  to  (l/MJ  , 
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hoory,  The  only  experimental  datf»  on  tie  variation  of  'A h-Si  with  ii 


^responding  to  a p of  five  certainly  does  not  fit  the  experimental 


polyvinyl  chloride  increases  Hmtx  from  0.10  to  about  0,17  and  radically 
decreases  the  width  of  the  dispersion.  The  value  of  p necessary  to 


than  five.  It  therefore  appears  that  the  predicted  frequency  of  HmaX 


six  but  more  than  one.  There  is  some  indication  from  the  experiments 
of  Fuoss  that  this  may  be  the  case.  Until  data  are  available  for  an 
unplasticized  polymer , however,  it  is  necessary  to  conclude  that  the 
theoretical  molecular  weight  dependence  for  a p of  five  does  not  agr'.o 
with  the  available  experimental  results. 


time  dependence  of  the  mechanical  properties  of  polyisobutylene  nv/" 
also  be  used  in  the  diffusion  theory  of  dielectric-dispersion,  '>*, 

corrected  dispersion  agrees  in  its  frequency  and  temperature  clrt.f..?.r‘ence 

' ■ t '■  F; 

with  experimental  results  for  polyvinyl  chloride,  bui  does  n'-t'/f Aye  th©-.,. 
proper  dependence  upon  the  molecular  weight. 


16.  R.  K.  Fuoss,  _.T,  Chew.  Sor,  63^,  2401  (1941). 


pointed  ou4 , hwpvr-r . that  the  plasticization  of 


describe  the  frequency  dependence  of  H(u>)  would  be  quite  a bit  less 


for  the  plasticized  polymer  would  vary  as  a power 


less  than 


In  summary,  the  resistance  constants  which  give  the  prop" 
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VIII.  SUMMARY 

If  an  irrepressible  noncross-linked'polymer  is  extended  in 
s^ch  a way  that  all  parts  of  the  material  are  strained  equally,  the 
velocity  v of  the  molecules  relative  to  each  other  may  be  expressed 
in  terms  of  the  rate  of  extension  cU/'ct  (Equations  (4)  and  (5)). 

The  configurations  of  the  molecules  at  uny  instant  may  be  specified 
for  our  purpose  by  £ .,  the  probability  of  each  configuration.  If  the 
temperature  is  above  the  polymer's  glast  transformation,  the  bond  lengths 
and  bond  angles  are  constant  during  extension.  Then  the  only  variables 
on  which  ^ depends  are  the  two  directional  angles  of  the  center  bond 
of  the  molecule  and  the  azimuth  angles  along  the  chain. 

If  the  molecules  are  not  cross  -linked,  the  intermolecular 
forces  during  extension  are  viscous  in  nature.  The  force  on  a single 
element  (one  skeletal  atom  together  with  its  side  groups)  is  propor- 
tional to  the  velocity  of  that  element  relative  to  the  surrounding 
material.  In  othor  words,  the  force  is  equal  to  , where 

is  the  resistance  constant  of  the  i«th  element  along  the  chain; 

'ir.  is  the  velocity  of  the  element  relative  to  the  center  of  mass  of  the 
molecule.  By  substituting  these  forces  into  the  equation  of  continuity 
for  £ , the  diffusion  equation  is  obtained; 

V ■ D-  [ Vf  + f W/kT  1 = |f . (20) 

The  diffusion  tensor  D may  be  calculated  from  the  by  the  method 

of  Kirkwood  and  Fuoss  (Appendix  I).  If  rotation  of  the  bonds  around 
their  azimuth  angles  is  completely  unrestricted,  the  potential  V/  is 


’ wyrrf > f ’ffl,  Try***  "* 
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gnttc»l  + q 


V = 


4 


(14) 


where  /o.t  and  z;  are  the  cylindrical  coordinates  of  the  i-th 
element  relative  to  the  center  of  mass  of  the  molecule.  We  may  con- 
clude that  the  molecules  change  their  configurations  by  diffusion  under 
a potential  proportional  to  the  rate  of  extension  de/dt  . For  a 
given  extension,  the  diffusion  equation  can  be  solved  for  the  probability 
of  each  molecular  configuration. 

The  rate  at  which  energy  is  supplied  to  the  polymer  may  be 
calculated  from  the  potential  and  the  probability  of  each  configuration 
(Equation  (27)).  Since  this  energy  is  also  equal  to  the  product  of  the 
external  stress  and  the  rate  of  extension,  it  is  possible  to  calculate 
the  stress  required  for  a given  extension  of  the  polymer.  A linear 
relationship  between  the  stress  and  the  strain, in  terms  of  a distribu- 
tion of  relaxation  times  £('*')  , may  be  obtained  from  a * ingle  solution 

of  the  diffusion  equation.  The  general  method  for  finding  E(y)  is 
summarized  on  page  31, 

Without  solving  explicitly  for  the  distribution  of  relaxation 

times,  it  is  possible  to  derive  the  dependence  of  the  stress  s(t) 

upon  the  temperature.  For  relaxation  at  constant  extension,  a plot  of 

-A/RT 

the  calculated  s(t;/  ' versus  tie  is  independent  of  the 

temperature.  This  dependence  agrees  with  the  experimental  behavior 
of  high  polymers  such  as  poiyisobutylene. 

To  find  s(t)  explicitly,  it  is  necessary  to  determine  the 
diffusion  tensor  0 and  solve  for  L ('>')  , To  get  a first  approximation 
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to  D , it  may  be  assumed  that  the  resistance  constants  are  all 

equal.  This  approximation  is  equivalent  to  assuming  that  the  motion 
of  the  polymer  molecule  is  the  same  whether  or  not  the  surrounding 
material  is  polymerized.  The  sit)  calculated  from  these  decays 

too  rapidly  in  time  to  fit  the  experimental  stress  relaxation  of  poly- 
isobutylene (see  Figure  3).  Moreover,  the  predicted  dependence  of  s(t) 
upon  the  molecular  weight  fails  to  agree  with  the  experimental  results. 

These  disagreements  are  not  surprising  because  the  calculated 
E(r)  applies  only  to  a dilute  solution  of  a polymer.  In  order  to 
determine  the  mechanical  properties  of  an  unplasticized  polymer,  it  is 
necessary  to  take  into  account  the  chain-to-chain  entanglements.  We 
make  the  jul  hoc  assumption  that  may  be  set  equal  to  €0  A , 

where  A and  p are  constants.  When  p is  five,  the  shape  of  the 
theoretical  stress  relaxation  agrees  quite  well  with  the  relaxation 
of  polyisobutylene.  Over  nine  decades  of  time  the  maximum  error  is 
five  percent  (see  Figure  5).  Thus,  it  is  possible  to  find  a diffusion 
tensor  which  gives  the  correct  time  dependence  of  the  mecnanical  prop- 
erties of  the  polymer.  The  dependence  of  sit)  upon  the  molecular 
weight  is  not  correctly  described  even  after  the  resistance  constants 
have  been  changed . 

The  diffusion  tensor  which  gives  the  proper  stress  relaxation 
of  a polymer  jnc„  also  be  used  to  find  its  response  to  other  types  of 
forces.  The  only  restriction  upon  the  applicability  of  0 is  that  the 
force  must  i esult  primarily  in  rotations  of  the  azimuth  angles  of  the 
molecules.  One  force  of  this  type  is  an  electric  field  of  audio 
frequency,  Kirkwood  and  Fuoss  have  calculated  the  dielectric  dispersion 
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of  polar  polymers  by  techniques  from  which  the  above  methods  have  been 
derived.  However , they  use  the  diffusion  tensor  applicable  to  dilute 
solutions.  It  is  possible  to  carry  through  their  computations  with 
the  new  diffusion  tensor  descriptive  of  an  unplasticized  polymer.  The 
results  may  be  expressed  as  a distribution  of  (electrical)  relaxation 
times  given  by  Equation  (102). 

The  corresponding  dielectric  dispersion  is  plotted  in 
Figure  6 against  the  frequency,  together  with  published  experimental 
results  for  polyvinyl  chloride.  The  agreement  between  the  theory  and 
the  data  is  quite  satisfactory.  It  is  believed  that  the  agreement 
would  be  even  better  if  the  diffusion  tensor  were  derived  from  the 
stress  relaxation  of  polyvinyl  chloride  itself,  rather  than  polyiso- 
butylene. The  predicted  temperature  dependence  of  the  dispersion 
also  agrees  with  observation.  The  theoretical  variation  with  the 
molecular  weight  does  not  agree  with  what  experimental  data  are  avail- 
able in  the  literature. 

It  may  be  concluded  that  the  ad  hoc  3et  of  resistance  con- 
stants correctly  describes  the  time-dependent  properties  of  a noncross- 
linked  polymer  at  temperatures  above  its  glass  transformation,  except 
for  the  dependence  on  the  molecular  weight. 
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APPEND n I 


DERIVATION  OF  THE  DIFFUSION  TENSOR 


The  resistance  tensor  for  a polymer  molecule  has  been  obtained 
by  Kirkwood  and  Fuoss  [7],  Their  derivation  is  reproduced  here. 

The  Cartesian  coordinates  erected  at  the  center  of  mess  of 
the  molecule  are  the  unit  vectors  e,  , et  and  e,  , with  ©3  in  the 
direction  of  extension.  The  elements  of  the  chain  are  numbered  from 
I to  2n  ; denotes  the  unit  vector  in  the  direction  of  the  bond 

from  element  k-i  to  element  k,  If  the  cnain  is  very  long,  it  may  be 
assumed  that  the  center  of  mass  of  the  chain  is  always  at  its  average 
position,  the  midpoint  of  the  center  bond.  Then  the  vector  , which 
denotes  the  position  of  element  k relative  to  the  center  of  mass,  is 
equal  to 


Xk 


k > 

k * * 


(103) 


where  a is  the  bond  length.  The  angle  9 is  the  angle  between  ©3  and 
, while  4 is  the  angle  between  the  plane  containing  and 

a-hrl  and  the  plane  containing  e,  and  This  definition  of  may 

be  written  as 


4>  = ?,)(?,.  §,)• 


(104) 
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In  the  same  notation,  the  azimuth  angles  X to  X.  are  defined  by 


X = 

m 1 

(9hw  » X?n»i  > 

(105) 

III 

c 

7< 

(2-*  > 

(106) 

f ( S-k-» i » 9-w  )( 9-k  * 9“K-i  ) 

2h-i  ?k  > h + Z 

III 

JC 

7^ 

1 (ak1  S-v..)(^,»  2w)- 

n -i  ? kj  2 (107 ) 

The  supplement  of  the  bond  angle,  q , is  the  angle  between  any  pair 
of  consecutive  vectors  9-w  and  9-k.,  • In  carbon  chains,  is  approx- 
imately 71°.  For  the  purpose  of  simplifying  the  following  calculations, 
however,  ri  will  be  set  equal  to  90°,  It  is  believed  that  no  serious 
errors  are  introduced  by  this  difference. 


The  velocity  of  the  element  k relative  to  the  center  of 

gravity  is 


k-i 

* Y 0/  hh-r,)» 

(108) 

jrh*l 

where 

O =.  a i 
~j  '3  J 

(109) 

and 

fl  — -(sin^)  0 e,  + (cos  <f>)  0e,  + (sin 0)  <£>  e^ • 
**  o 

(no) 

If  the  molecule  is  moved  through  a quiescent  fluid,  the  force  on  the 
element  U is  £k  yk  , and  the  torque  at  the  position  of  another  element 
L is  Ck  (fk- JTi)  * Yu  • The  total  vector  torque  acting  at  the 
position  of  element  i due  to  the  forces  on  all  the  elements  is 
Zn  i 

T = Y.  x Tw  ~Y_  vk-  ^ (lli) 

k - 1*1  U*  *>♦! 

The  component  of  X in  the  direction  of  9-;  is  the  scalar  torque 
Tt  associated  with  the  angle  X^  . Substituting  Equation  (108)  into 
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Equation  (ill)  and  rearranging  after  dotting  with  a-  , we  have,  for 


Jn-I 


2n 


* *r 
Z n 


where 


and 


k » i * I 


} aTj  +/°u-Q-o  > 

(112) 

if,  OL.  • Y;  k'  Q-  • 
~ i ~ lj  * i 

1 

-Mu  -*j  - y 

h+l  4 

i<Ul 

*•  w»j«i 

iSn 

(113) 

S-I-Xiow  ~^T  Cw^*?L.k. 

(114) 

krn+i 


The  dyadics  are  defined  as 

s (rK-rj xVr.)I  ~ (v^Urv-rO^  (115) 

where  I.j  = . The  dyadic  denotes  with  identically 

zero.  A corresponding  result  may  be  obtained  for  the  lower  half  of 
the  chain.  The  torque  associated  with  the  angles  6 and  ^ is  the 
total  torque  on  the  molecule  referred  to  the  center  of  mass  : 

To  ~ £oa'Q°  + Y ' (116) 


where 


2n 


e 


k=  j+> 

j-> 
r~ 


A f * -a-. 
7 ^ j 

I 


M -f  I 


; <.  n 


(117) 


/?«  “ T ^ (-k'-k?" -k-K  )* 


k=« 


and 


(118) 


I 


f 
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The  dyadic  denotes  with  f t identically  zero.  All  the 

terms  of  the  resistance  tensor  p are  given  by  Equations  (113),  (114), 
(117),  and  (118). 

It  is  prohibitively  difficult  to  find  the  inverse  of  p and 
thereby  to  obtain  0 = . However,  the  diffusion  tensor  may 

be  approximated  by  kaT(/5)  , where  p is  the  average  of  the  resistance 

tensor  over  all  possible  configurations  of  the  molecule.  Substituting 
Equation  (103)  into  the  definition  of  k , we  may  show  that 

a-1  [ I - (a-.*  1 ? i + 1 


a.  ■*..  -a.  = 

~ i - » jK  - j 


V' 


0. 


The  value  of  (9v  is  1,  0,  1/2,  1/4,  3/8,  5/16,  ...  for  equal  to 

l,  Ut,  1+2, . . . , and  approaches  1/3  as  becomes  large.  Without 

T ^ 

much  error  in  p , (at*a  ) may  be  set  equal  to  1/3  for  all  Using 

Equation  (113).  we  then  have 

* O ^ 


p-r 


f 4aS;I  CkC'-o 


i 


s n 


(120) 


The  value  of  is  given  by 


OO 

2 h 


(t 00  — -y  2_  ^ ~k’  ±-  • 

ki  1 

Using  Equation  (103),  and  neglecting  the  term  -i-  9-^^,  , we  have 


(121) 


2n 

k-* 


vrK  = 2a 


z v- 


/ V ^ <?.  a • q 

Z_  2_  Snn  ' 

K*  1 k-M  >r jr  k*i 


so  long  as 


<0=  c 


2h*  1 - k 


(122) 


(123) 
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From  the  faot  that  ^ > 


k=i 


Also,  it  may  be  shown  that 

O-  - p . - o. 

iZ  *3  f " Oi 


(124) 


(125) 


The  components  of  the  diagonal  diffusion  tensor  are  readily 


found  from  the  average  resistance  tensor,  and  are 


Pee  = DohV4^h) 

(126) 

and 

0 — D n &.  /?  q( 0 , 

*‘j  3 

(127) 

where 

0o  = 3k.r/^n\ 

(128) 

and 

(j~ 

Ws  !♦  1 

i ? n-»  1 

i.  S’  O 

(129) 

If  it  is  assumed  that  !^k  is  equal  to  t ^(O  is  approx*- 
imately  £n~  /Z,  The  nonzero  terms  of  the  diffusion  tensor 


are  then 


and 


%*=  o./2 

p;l  = 


(130) 

(131) 


The  diffusion  operator  is  equal  to 


v-ov  = &)♦  *7? It*  I (02) 
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and  the  eigenvalues  and  normalized  eigenfunctions  of  this  operator  are 

2 n-l  _•»  •> 

‘•mi  4 2 


x=  naM-iHn-iW 

^ i=2 


where  P(  (cos  ©)  are  the  associated  Legendre  polynomials. 


(133) 


(134) 
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EVALUAT  ION  OF  dxdx 


From  Equation  (70),  dx  is  equal  to 


when  £•  - . Using  Equation  (103)  of  Appendix  I and  neglecting 


Z ~k*>  > 


vg  find  that 

2 n 


2n  J 


is) 

2h 


I r K4-2  W -ht2. 
2 it  ; 


(135) 


Yi  Y1  ~i'~k  + Za*  H y~  (2*+t- j)gj‘Su  (136) 

j k*, 


and  Y~  Z?  Y~  (j-i)(g^9iXg,-?J-t-2o.2  y~  y~(2ntl-jXeya  )(e3-ak).(137) 

i'.:l  jr^  W-j’i  jin»2  krn»i 


From  r2-  3*‘  , it  follows  that 


X>')(£,+  £1  !>*'-;}  p/;\ 


i-2  I' - i 


jrh«2  K-nt-i 


where 


C=  ^ f [v  *<■'  3<?J-9j>(S3-  M £d<t 


(138) 


(139) 


D(x> 

In  order  to  find  r -k  , it  is  necessary  to  express  O-j  in 
terms  of  the  angles  0^0,  X2 , From  the  geometry  of  the 
molecular  chain,  it  may  be  determined  that 

~J’  1 A(ho^h-,)' •••■  ~j  jih  *140* 
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where 


A(Q,^/  = 


/ cos  0 Cos  <fi  -sinj>  SinQcoi^ 

j CoS  Q sin  coS^  sif>6sm^ 


- sin© 


(141) 


It  will  be  assumed  that  ^ , the  supplement  of  the  skeletal  bond  angle, 
is  equal  to  90°. 

For  j?  n+  I , substitution  of  Equation  (140)  into  Equation 
(139)  gives 

n(X^  r~f tt(  <°>  <***•>  (h*l>  <>*-•)  ^w>  s{k,,)  .A  i 

: /*■  ; V : v>  ^ :Kj^  (142) 

where  K is  the  dyadic  fse  , and  the  fourth-order  tensors  yu°, 
and  are  defined  by 

v(e<^s * (143) 

= ’ (144) 


and 

nc<sr&  = A ? yS  • 

(145) 

The  scalar 

function  U of  any  dyadic  JJ  denotes 

u (T)  * T„ + F«  * 2 IT3J  • 

(146) 

The  double 

dot  product  is  defined  as 

lAS  J^r  ^ £««v«  ^ep<S  . 

c„  *» 

1 

(117) 

and 

(148) 

Substituting  for  ^ from  Equation  (134)  of  Appendix  I,  we  obtain 


p(A)  =U(n"!):  |-T°:  N(k\N  N^*0 : K 

i*1  n 3n.| 


x tt  s:  n c , 

vx  i*j 


(149) 
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tr  2xr 


where 


U(MW>  = JLj  J Ut^Y^^^^neded^, 

©•O 


iTT 


M1''  = — f J* 


2 IT  ] ^ 

o 

j?ir 


(150) 


(151) 


(152) 


[y  = A^i-f + r *«»? 


(153) 


»/2 


and 


Ym‘(e  4)  = P'Vco*©' 

V ' 4ir(JM«0)!  J r* 1 ,e 


As  the  first  step  toward  finding  R*  * , it  is  possible  to 

evaluate  the  product 


(154) 


L « n“":  n""":...:  n“"’:K. 


(155) 


Integration  over  gives 


N 


<V 


- [ Ao  + f + x A_  j I, 


(156) 


where  Aq  = - e}e(  and 


A*  = ?,  ?.  + ?,§3  * s2  ?3  * ?,  ?,  /£T 


(157) 


By  substitution  into  Equation  (155)  and  contraction  of  the  dot  products, 
we  find  that* 


L =(-if""  (i)'(f)'..  n [ f>  - (*v  ,(«),<•>;.,) 


(158) 


^Equations  (158),  (164),  (165),  (171),  and  (233)  may  be  proved  by  induction. 
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where  (°)^  , <0,^  and  (0^  denote  , (&m  + and  )//7  , 

respectively.  Later  in  this  appendix,  U)  and  C2^  will  also  appear; 

they  denote  (S*  + Sm  ) and  (S(n  - )//*" i , respectively.  The 

1 ^ % \ 

quantity  <r  , defined  as 

ffS5I(0v  (159) 

k 

is  the  number  of  m,  that  are  ii  for  . L is  zero  unless  ail 

the  are  O or  + ! . m = **  » no  consecutive  m.  ?nd  m^+|  Kre 

zero  simultaneously.  If  L is  not  zero,  its  absolute  value  depends  only 
upon  <r  . 

The  next  step  in  evaluating  P.^  is  to  find  PI1'5  and  its 
double  dot  products.  From  Equations  (144)  and  (151),  we  may  show  that 
M**1  is  the  sum  of  and  , where  is  symmetric  under 

transposition  of  the  first  and  third  indices  and  under  transposition 
of  the  second  and  fourth  indices,  whereas  n.  is  antisymmetric  under 
either  transposition.  From  the  nature  of  the  IT  function,  it  is  pos» 
sible  to  show  that  only  the  symmetric  parts  of  and  L contribute 

to  P™  : 

JK 


P*A)  =U^MC*’) 

1 jk  - - - s -* 


h*0  ..Owi' 


. I 

Vs  * ts 


% 


* n i 

n 

Vi 


(°V 


The  value  of  obtained  from  integration  of  Equation  (151)  is 

= [A  a„4a,a.  HA-A.]S^|(a.a.^a.aJS^ 

Rearrangement  of  the  terms  and  subtraction  of  the  antisymmetr ic  part 


(160) 


(161) 


I 


r 


APPEND  IX  II 


• Kxl  iJt 

give  H,  , which  may  be  expressed  as 

[n;vL,s  = + (&-,  Pp  * -E..rps-n'“)v 
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(162) 


/C2^  (0^  /O  0)  -0)  \ 

An’  = y(  e>;  -*>,  m;i,  ?<v  - f ( o»;  , 

H k W -n  ev 


/ 0 (.,)'  *0) 


where 


(163) 


is  zero 


C =(§,§,-»- et?a)/2  , 0=  , E = ejel,andF=eei. 

unless  =0,  ±1,  or  ±2.  Also,  it  may  be  seen  immediately  that  Mj'*  * 

is  zero  if  either  m^=  O and  * ±\  , or  vice  versa.  Hence,  the  total 

..(h-H)  -.(k-O 

product  Hj  rlf  (except  perhaps  near  its  ends)  may  be 

divided  into  a set  of  scalar  products  of  the  form 

B‘V:  A<V"’>  , .fere  |m,|  = | | = 2 , 

and  all  m,  in  between  (there  being  v of  them)  have  the  same  absolute 
value  0,  1,  or  2.  Each  of  these  scalar  products  may  be  contracted 
separately,  to  give: 


/-II 


V2 

V*- 


<2>v— 

mv« 

=0,r>l 

(0vM~. 

|«v 

\ = \tr>  : 

• 

i>l 

\=ls>°  t/,ri4) 

The  connection  between  the  M,  and  L,  can  be  expressed  as 


the  product  of  L,  with  the  f1s  back  to  B 


( It-*--  i) 


, where  k-w-l  is 


*Note  that  the  indices  in  the  right  side  of  Equation  (162)  have  been 
transposed  from  their  order  in  Equation  (144). 


I 
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the  largest  integer  such  that  |rnk_r_(|=2  . The  nonzero  values  of 
this  product  are 


8“ : M?— >, 


:L,= 


(°V  C'C 


- y(-  yffi.,  ('C  H =l,^o 


Where  W,  from  Equation  (158),  is 


w-(-oi'k""(if(-l1  FT  [«\- I'LioV’vO- 

V k*' 


From  the  definition  of  1 in  Equation  (143),  it  follows 


that 


y(/i“>=  t 


coS2e--j  o 


sin28  \ 


o 

sin  28 


O -cos2©  - -j 


( 


Substituting  U(pr ” ) into  Equation  (150),  we  find  that 

V(vM>  - ^mc  [ K £ * ~ + t £1°*  htrr  ^ 5 ] » 

odd  a 

where  G = - ee  -t- e,et } H=  e3e,  , J=  e, e,  + €>,  e ( - 

The  numbers  a v are  the  expansiun  coefficients  of  sin20  in  terms  of 

the  Legendre  polynomials  of  cos  © [17]: 

sin  20  = (.Cos ©J ^ 

odd  a 


17,  W.  E.  Byerly,  "Fourier  Series  and  Spherical.  Cylindrical,  and 
Ellipsoidal  Harmonics,"  Ginn,  1893,  p,  183. 


(166) 


(167) 


(168) 


(169) 


i 
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with  a(  = 3tt/8,  and 


_ (w-2X*-+2)(2«+5) 

~ ^*l)(«  + SX2f  + I)  a*' 

The  connection  between  UCM^)  and  the  can  be  expressed  as  the 
product  of  UCM^)  with  the  M,  up  to  /\(-r'*r~f''  * where  n+r+l  is 

the  smallest  integer  such  that  |mntr<(|=  2.  . The  only  values  of  this 
product  which  are  not  zero  are 


(170) 


+ r*0_ 

- (2)*„  s;. 

A~O,r-0 

A(-i)"<0L,<0)..,-<^. <2) s~. 

1-2, r>l 

£-2 , Y — 0 

- (-t/S  ISr  ('KM.,- CL,J'L 

» » rt  .1 

odd  Jt,  r * 1 

(171) 


The  above  connection  formula  completes  the  determination  of 
P^X'  . For  a given  , R^  may  be  evaluated  by  substituting  H, 


m0  , and  the  tr>  into  Equations  (164),  (166),  (166),  and  (171).  It  may 

0<>4  . , 

be  seen  from  these  equations  that  is  zero  unless,  for  j > k ? t-> -*-l , 


1) 

2) 

3) 

4) 

5) 
S> 

7) 

8) 

6) 
10) 


m0  = 0 

^%-0  for  a*n 

+ ^(2)h+i  = I 

[m^|=0j,ot-2  fo  r n-l«jc^k-! 

C'\(0)v,+  C°\(')ri  ~ O fo*  n+t*  k-2 
(<>)„.,  (o)k  + (0  J--\  i-  (*W<0k  - C2)u.,0)k  = I 

|mj  = 0 o>-  ! for  j-  2 

(°)aCo)r=  0 ^or  -3 

K-H 


APPENDIX  II 


71 


Corresponding  selection  rules  occur  when  jiki-n. 

After  Pj^  has  Deen  determined,  dx  may  be  found  from 

Equation  (138).  From  selection  rules  9)  and  10)  above,  it  is  apparent 

that  there  is  at  most  one  nonzero  term  in  the  sum  over  j , corresponding 

to  the  largest  j for  which  jrrv|  = | . The  other  selection  rules  permit 

more  than  one  nonzero  term  in  the  sum  over  k only  if  there  exists  a n 

such  that  jm  | =|rv>^,|  =...=  Im^.^1  = I,  | > | = 0 , and  <\  is  the 

largest  integer  for  which  |n\J=2  . For  this  case,  there  are  r nonzero 

terms  in  d*  , It  may  be  shown  that  the  cross  product  terms  in 

( with  i<  7s  k ) contribute  nothing  to  -Q(X).  Furthermore, 

D(>0 

all  of  the  squared  terms  rjX  are  eoual  under  these  circum- 

stances, so  that  we  obtain  the  correct  £"1(\)  if  we  use 


d>d;  = 4«v»>-j)’R  rc 


(172) 


The  factor  w is  present  only  if  there  exists  a satisfying  the 
criterion  above.  Since  r changes  £1  (X)  very  slightly,  it  will  be 
neglected  in  the  later  calculations. 

From  the  equations  for  , it  may  be  seen  that  dK<dA 


does  not  depend  explicitly  upon  the  exact  values  of  the  nonzero  , 
but  rather  upon  certain  functions  describing  their  number  and  relative 
positions.  In  order  to  express  c^A dA  succinctly,  let  us  denote  by 
fx  the  number  of  rn  between  and  that  are  equal  to  zero: 


h 


k — \ 


(173) 
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Similarly,  0 and  C denote  the  number  of  ± l's  and±2's  respectively: 


k-i 


k-i 


d=Y_  I and  $ = Y1  tZ\‘  (174) 

<^=n+l 

Also,  let  « , (i  , and  y denote  respectively  the  number  of  sets  of 
0's,  ±l*s,  and  ±2’s  between  »■«„,,  and  mk.,  , where  a set  of  0's,  for 
example,  is  an  uninterrupted  string  of  one  or  more  rn^  having  the  value 
zero: 

Ui 


« s Y.  [<°\-  (oV0,v] + (°v, 


(175) 


?=£  K-om,,] 


k-t 


K’Vrl- 

V"*' 


(176) 


(177) 


Finally,  let  cr  be  defined  by  Equation  (159)  above.  Then  d^d*  depends 
only  upon  j > and  upon 

and  | mk  | . As  an  example,  let  us  choose  M.~  ^ k_,~  O . If  dxdx 

is  not  zero,  it  is  necessary  that  |m01J:=  2.  , mk  = o , and  jfn(Ul|=l. 

Also,  fn0=o  , and  , For  this  case.  Equations  (164),  (165), 

(166),  (171),  and  (172)  give 


d d 

X A 


- Ik. 

- j**'' 


2-2*  2T  , Su  , . ii  .2?  . 2<T 

(i»)  (4)  (j)  , 


(178) 


where  j is  the  largest  integer  such  that  m.  ^ 0 . Values  of  £, 
lmk-i  | , and  1 1X1  u | other  than  those  chosen  above  give  the  same  dependence 

of  dxdx  upon  , 0 , , /a , 0,  ^ , & , and  j,  but  with  a multiplicative 
constant  different  from  16/9. 
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EVALUATION  OF  fKX) 

From  Equation  {12 ) , -O  (^ ) is  approximately 

N(X)dj<,K'i  (179) 

x C-A« 

where  N(x)  is  the  number  of  havirg  some  particular  absolute  value, 
and  X is  the  average  eigenvalue  of  these  d,  . From  Appendix  II,  the 
''absolute  value  of  dx  depends  only  upon  , |mh+il,  |n\-,K  and  l^kU 
and  upon  *,?/)/*)  ^ (j  ffj  and  j . Hence,  N(X)  is  the  total  number 
of  nonzero  dx  that  may  be  obtained  by  rearrangement  of  the  while 
keeping  J?,  |mn<(  . j constant.  The  calculation  is  simplified 

somewhat  by  introducing  the  variable  s in  place  of  j ; s is  equal  to 
the  total  number  of  m from  up  to  nrv_(  } inclusive: 

S=  (180) 

The  dependent  variable  w will  be  used  to  denote  ^L  + ^+C. 

As  an  example  of  the  method  for  finding  OCX)  , let  us  assume 
that  J.  = fo(>_|  = 0 , and  thus  that  jmn<||=2  f my*  O , and 

Evaluation  of  X/(X)  may  be  facilitated  by  the  following  theorem: 
there  are 

possible  ways  of  placing  A indistinguishable  objects  into  & distinguish- 
able sets  with  at  least  one  object  in  each  set  L18).  Between  mn+1  and  rr>K 


18,  E.  Ising,  Ze.it s . £ Physik.  31,  253  (1925),  Footnote  1, 
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there  are  thus 


(KXmXE) 


arrangements  of  the  into  their  sots.  The  factor  2^**  is  the 
number  of  possible  arrangements  of  the  plus  and  minus  signs  of  the 
It  is  also  possible  to  show  that  there  are 

(I:!) 

different  ways  of  arranging  the  order  of  the  sets,  keeping  |mh4||  = 2, 
mh.,  = 0 , and  with  alternate  sets  formed  of  ±2*s.  Finally,  there  are 

2r 


arr 


angements  of  the  ±l's  and  0*s  between  and  , keeping  (0)  (o)+)=0, 

■ i 

. Combining  these  results,  v?  have 


mL*0  . and  I m.  . | = 

••  » • j-i  • 


(181) 


The  additional  no'ror  two  appears  when  both  sides  of  the  molecular 
chain  are  included  in  N(\)* 

From  Equation  (133)  in  Appendix  I,  the  value  of  X is 


2n-| 


x = os  Y_  [•  - IV1!] 


when  1-0.  Changing  the  summation  to  an  integral: 


W W4| 

, <r  f JCv^)  1 
-L  w J + s J (,_  A-,*  I 


v 


or 


X = 


oD.fo-n  + s)J 

(«- w)(r.-w  -?) 


(182) 


(183) 


(184) 
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From  Equation  (178)  in  Appendix  II  and  the  above  results,  it  follows 
that  the  contribution  to  jQ(X')  for  j£=mk_(  = 0 is 

2p.«  r-2T 


= H ¥• (+f  (4-f 


<r,s,u.,0 

f r 


(185) 


where  When  the  sura  is  taken  in  the  order  <r , s , 0 ^ ? , “ , 

and  IT  , the  limits  are  approximately  s/2  4 f ^s,  s tn-w, 

(3*0$  n-t-*,  U h-r,  C*«*1T,  0$Y$n/2, 

and  <Tn  -h  (4-  (j  0)(r*- w- »)  $.  ^n-w)(n-w»)/riPs. 

If  the  molecule  is  very  long  (n»l  ),  the  summation  may  be 
evaluated  quite  accurately  by  replacing  the  binomial  coefficients  with 
their  asymptotic  expansions  for  large  numbers  [19],  As  an  example, 
the  first  sum  is 


<r-  l 
s-cr- 


1 


). 


<r 

with  limits  s/2£rr£s  and  cr  $ (*-w  -s)[[X(n-w)/tiD,j  - (4$+0)}/n. 
By  expansion  of  the  binomial  terra  for  o'  near  25/3  , the  sum  may  be 
replaced  by 


(2s-3<r)2 

2c 


This  integral  is  approximately  1/3,  provided  that 


s $ 


3(*-~)[hVn~w)‘(4*  + ^l 
2ntf  + At^n'w)-C4?  + 0)]}  ” 


(186) 


19.  K,  Jeffreys,  "Theory  of  Probability,"  Oxford,  Oxford  University 
Press,  1939,  p.  50. 
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If  s is  greater  than  this  quantity,  the  integral  is  approximately  zero. 

In  this  way.  Equation  (185)  may  be  evaluated  by  a sevenfold  integration. 

The  result  of  the  integration  is  the  partial  contribution  to  fKX) 
for 

[nwl.  = 

where  C s 3X/2nD0.  (188) 

To  obtain  the  total  value  of  XI  (.X)  , it  i3  necessary  to  find 
X , N(X)  , and  d-d!  for  all  values  of  and 

|mw|  which  give  nonzero  *Jxdx  . Change  of  these  variables  alters 

by  a multiplicative  constant  (see  Appendix  II),  alters  (3  slightly 
in  terms  of  V and  «,  and  affects  the  term  ) in  N(\)  . Further- 
more, the  value  of  X is  dependent  upon  Jl  , but  this  effect  can  be 
neglected.  When  all  the  partial  contributions  to  -0-(.A)  are  totalled, 
the  final  result  is  thirty  times  the  partial  sum  given  above; 

n(x)=  30[n(X)]x-o  . dag) 

*v, 50 

Frc2  lGi*3  \ 187)  and  (189),  we  can  show  that 

_ . . A 4 4 

XI (*)  — n o-  • 

JX(oo)  may  also  be  obtained  from  Equation  (48): 

XK“)  = [ -2e?)]2, 

i 

that  is,  H(oo)  is  the  average  over  all  possible  configurations  of  the 
square  of  the  sum  over  L of  (p’-2»;)  , By  using  Equations  (136)  and 

«137),  it  may  be  shown  that  Equation  (191)  leads  to  the  same  value  for 
O(oo)  as  is  given  above.  This  agreement  serves  as  an  excellent  check 
or,  the  magnitude  of  O(X). 


(190) 


U*j.) 


- C- \U  -4Sa(crl)  + 4-  577^,}  , 


(187) 
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Fro n Equation  (SO),  the  force  at  constant  extension  is 


s(t)|  = cj 

* o 

If  the  function  Z(c)  is  defined 


e ^ E.  (*r)  eh- . 


Z (c)  = — — f < J42  -+-  3c  + 24|rv(c+l)  - ^c*ll)  - *£2S*}2}  _ §i£l3)  ) 

1 ' (c»3)  l 1 (c*<)  Cc-M)*  ic+l)3J  ' 


Squat  ion  (75)  may  he  written 


as 


(192) 


(193) 


t(ri=  <“«> 

“a  S* 

rhere  c~  \/r  and  T*  = n ?.  / 2 IqT  . Substitution  of  into 

Equation  (192)  gives 


s(t)j  =4N.,V.TcfzCc)e*(t/T-Kd(^)- 


(195) 


he  integrand  has  teen  evaluated  numeric  ally  as  a function  of  l/c 
ntegrated  graphically.  Ihe  resu.1+‘I 


and 


Table  I 

Theoretical  Stress  Relaxation  for 
t/% 

?£=  C. 

s(t)/sCo) 

0 

1.753 

1.000 

0.0035 

1.676 

0.956 

0.01 

1.553 

0.886 

0,035 

1.257 

0.717 

C;1 

0.846 

0.483 

C.35 

0.310 

0.177 

l.C 

0.G7C 

0^040 
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The  normalized  stress  s(t)/i(o)  is  obtained  by  dividing  the  integral 
ty  1.753. 

The  graphical  integration  nay  be  checked  by  finding  the  exact 
value  of  the  integral  when  r =o  - By  ordinary  analytic  nethods,  it  nay 
be  shown  that 

C ZC^dC  _ _7 
J c:  4 ’ 

O 

agreeing  with  the  value  obtained  by  graphical  integration  to  within 
0,2%.  Ihe  stress  at  t = 0 is  thus 

s(t>/  s(o)  is  plotted  versus  the  logarithm  ur  t/%  in  Figure  3, 

2 04 

with  "ro  equal  to  10  * hours. 


(196) 


(197) 
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EVALCATICK  OF  Q(h)  I CSL  gt=5.A 


— # i. 


It  is  the  purpose  of  this  appendix  to  calculate  the  diffusion 
tensor,  the  expansion  coefficients  cf  the  potential,  and  G(X)  that 
result  from  the  assumption; 


u 


C.A 


A(2«  + I-0P. 


ijn 

i ?>»*• 


(19S) 


The  derivation  of  the  average  diffusion  tensor  is  carried  out 
in  Appendix  I in  a way  somewhat  mere  general  than  was  first  given  by- 
Kirkwood  and  Fuoss,  and  its  results  are  directly  applicable  here.  Frcm 
Equations  (126),  (127),  (128).  and  (129).  the  only  nonzero  terms  in  0 
are 

— ss  _ £p*»Kf 
fee"  4A*f 


and 


where,  as  before. 


,A  Jn^TlI^nP1  ' 

— * n"  i m ♦ J 


(p*'Xp*2)P. 


(199) 


(200) 


°o  = 3VjfZm*<n 


(201) 


The  normalized  eigenfunctions  remain  the  same  as  in  Equation  (134),  but 
the  eigenvalues  are  row 


(p«»](p*2)Q.  r,  v r 'r * j 

K=  2«:r  ij  ->) 

The  expansion  coefficients  of  the  potential  are 


2e-P 


(?*2l 


(202) 


(203) 
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Using  Equation  (103)  of  Appendix  I and  neglecting  la*.*,  , we  find  that 


In 


n m 


2r*  j! 


~f^r<  - 2c-s  y + 2°i  y_  Y_ b(,J  -/-* 


(204) 


l z pk  Z Ik  » ** 


and 


j=2  fc.. 


Ji-  • 


where 


«)  = 


j ^ « 


(205) 


1- 


(206) 


14- 


V 

k*r‘ 

By  evaluating  b(j)  2nd  using  p*-  2a*  = r2-3z3  we  find  that 

2*  i 


r>  f> 


<k=^[z  si' •>•/■■ pf  + 1;  t: (2-'-r' o <*»> 


i* 2 v-r  j 


n>-»2  k.tn»J 


r is  unaltered  by  change  of  the  diffusion  coefficients  and  may  fc*a 

•i  ** 

found  from  the  formulas  in  Appendix  II.  For  the  example  given  previously. 
( SL  — rr» k ||=  O)  , substitution  of  r;fc  onto  Equation  (207 ) gives,  for  j ^ 

— I<  2'S  It  I.  2)  2f  }<r 


<U  * = ,fea:A‘  - 


‘MV*  A *»  X • c i 75 

« - {$>  (T)  M (|)'  {*>  & W- 


(20S) 


li  is  now  possible  to  find  Q (A)  . N(X)  has  the  sane  value 

as  in  Equation  (151),  but  X.  is 

2*-l 

-Ho-!} 


— (p*')(p+2)q,  \ f,  | — 

A ~ ~T^  — L "h 


(209) 


T 2 


when  A — O . Using  an  integral  as  an  approximation  tc  this  summation, 
we  have 


T_  J i 1 rr 

A 2 A 1‘  - i nn*i  * sl'i 


» 


a ^(h-vr-s)***  J P^11  J* 


j 

5, 
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From  Equations  (179),  (181),  and  (208),  the  contribution  to  -fl(X) 
when  A ■*  0 is 


2 . 


[o<%0=  L J 

(TS.ii.V, 

x(c-.:)(?::xfc'J(:::  )(»-,)• 


^ 3 f ♦ 2 r 


'V.*° 


(211) 


The  limits  on  this  sum  are  the  same  as  are  given  in  Appendix  III,  except 
that  the  last  one  (derived  from  X fc  \ ) is  now 


/ n-w  >P**  < t(n  f 2AX  _ I A f + 0 \ f I I 1 1 

(n-w-J  ' I + <r  \.(,p*aVD#  ' W 'L(n.w )f*'  r, P-* * Jj- 


(212) 


By  asymptotic  expansion  of  the  binomial  terms  for  large  values  of  the 
variables,  the  sum  may  be  tr ans forced  into  a sevenfold  integral  and 
evaluated  by  the  procedures  exemplified  in  Appendix  III.  Six  integra- 
tions yield 


"V-1*0 


g A 8«Va^ 

4-5  (?+>)* (2 r 3) 


i(c) 


(213) 


-/IfjLf)  c (P*3>  /2p-»3 1 

where  ti(c)  = Cc+3)  p*'  j (*-2)  771 ; [l  - x ( ]dx,  (214) 

I 


with 


c = 


(p*2) 


r;x. 


(215) 


"he  seventh  integral  is  not  generally  expressible  in  closed  form. 

As  before,  the  total  value  of  O(X)  is  the  sum  of  the 
contributions  for  various  j2.,  |nr,n,(|,  l^u-t  K and  |mw|  , The  rela- 
tive magnitudes  of  the  contributions  are  unaltered  by  change  of  the 
resistance  constants,  and  Q(X)  is  thirty  times  the  partial  sum  above. 
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Ft on  Equations  (213)  and  (214)  ve  may  shew  that 


Q{«0  - 


8AV 

3(p-*-02(p'r2){2p+3) 


(216) 


As  in  Appendix  III,  it  is  also  possible  to  find  Clip*)  fr&a  Equation  (48); 

a (e»)  = [ Y1  T f\  - 2 ** ) f - (217) 

t 

By  using  Equations  (136)  and  (137)  we  find  that  this  average  leads  to 
the  sane  value  for  as  is  obtained  above,  indicating  that  the 

magnitude  of  DO)  is  correct . 


S3 
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CALCULATION  OF  THE  CORRECTED  STRESS  RELAXATION 


When  p is  equal  to  five,  $(c)  in  Equation  (214)  becomes 

T 4 J2 

«(c)  =(c-3)~  i)  (<■*■  2)3  { > - * )d*. 


/OlO\ 


We  could  find  a solution  in  closed  form  for  this  integral,  but  it  would 
to  quite  involved.  Instead,  ve  shall  use  the  approximation  that 


**C  4 J3 


f(,~z)3*~fcJ.  S U-0‘f4-j^T-(|^]-2.4 40.  (219) 

Fron  Equation  (91),  £(7*)  is  then  equal  to 

e'W  = {7JS7  * T<‘-||nr  g&>  - f^rPi}.  (220) 

From  Equation  (94),  the  stress  at  constant  extension  is 

,(t)|(=  9.0018 17  X.- (.T « f ^7. 287 + (c. <)‘[2  (4-ig? - 3 ) -(gf  )Tj  k <221 ) 

o 

1 his  integral  has  been  evaluated  graphically.  The  amount  of  graphical 
integration  nay  be  reduced  considerably  try  analytic  app roxinat ions  to 
the  tail  of  the  integral  from  Ca  to  as  , where  cs  is  any  number  much 
greater  then  unity.  Also,  when  t/rj  is  nuch  smaller  than  unity, 

S(t)-  s(0)  nay  be  determined  analytically.  Values  of  the  integral 
for  various  t/%  are  given  in  the  table  below.. 


[ 
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The  graphical  integration  nay  be  checked  by  conput mg  the 
exact  value  of  5(0)  . For  arbitrary  p, 

fp-2)?  NI.nls.Tt  fi  J £(c) 


5{0)|  = 

U 3(9,\) 


3{p-‘)J(2r3)  l dr 


(222) 


By  ordinary  net hods,  ve  find  that 

f 3 c ft(c!  . _ (p*1/  (2  ?*2’' 

jC  <it  C _ 2(p^2Kp^3)  ' 


(223) 


For  p=5  , the  integral  is  53,79,  agreeing  with  the  value  obtained  ty 
graphical  integration  to  within  0.2%.  It  follows  that 


fc{pi»3){Z»*3) 


N.nk.Tc. 


For  $>  = 5^  s(t)/ 5 (0)  is  plotted  versus  the  logarithm  of  t/%1  in 

Figure  S,  with  v0  set  equal  to  10"* hours. 


(224) 
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Table  II 


Theoretical 

Stress  Relaxation  for 

Ci=  C.A  i.5 

t/Tj 

s(t)/s(o) 

0 

59.86 

1,000 

10“9 

57.05 

.953 

10"8 

55.71 

.931 

10"7 

53.72 

.897 

10"6 

50.88 

.850 

10-5 

46.51 

.777 

io"4 

39.56 

.661 

io"3 

29.40 

.491 

io"2 

15.73 

.263 

0.017 

12.62 

.211 

C.03 

9.30 

.155 

0.06 

5.77 

.0964 

0.1 

3.74 

.0625 

0.17 

2.14 

.0358 

0.3 

1.07 

.0179 

0.6 

Oe35 

.0059 

1.0 


0.12 


0021 
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APPENDIX  VII 

CALCULATION  OF  THE  CORRECTED  DIELECTRIC  DISPERSION 


The  general  solution  of  the  diffusion  equation  for  a polymer 
molecule  in  an  electric  field  was  given  originally  by  Kirkwood  and 
Fuoss  [7],  Their  solution  is  as  follows.  Let  be  the  unit  vector 

in  the  direction  of  the  electric  field.  Let  denote  the  expansion 

coefficients  of  e3  , where  ^ is  the  dipole  moment  of  the  molecule: 

and  Xj  as  before,  are  the  set  of  eigenfunctions  and  eigenvalues 
of  the  diffusion  opera-tor  V-Q-S7  , Let  O(X)  be  the  sum  of 
for  all  * , normalized  so  that  £!{«)  ~ l ■ 


I 


Z r*r‘- 


The  normalisation  constant  iff  ?3)z  is  the  average  over  all  possible 
configurations  cf  the  squa;  : cf  p-e3  • The  reduced  polarization  of 


the  polymer  is 


QM  -ftT^T  • 

J ( I + tco-rl 

O ' 


where 


x = i/r 


Kirkwood  and  Fuoss  calculate  Q(u>)  upon  the  assumption  that 
. Consequently,  their  result  is  applicable  only  to  a dilute 
solution  of  the  polymer.  It  is  the  purpose  of  this  appendix  to  solve 


(225) 

(226) 

(227) 

(228) 
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p 

for  Q(>»v  when  = . When  p is  five,  this  result  should  give 

the  polarization  of  an  unplasticized  polymer.  It  will  be  assumed  that 
the  polymer  is  polyvinyl  chloride. 

In  polyvinyl  chloride,  the  individual  dipole  moments  are 
associated  with  the  C-CH  bonds.  Let  be  the  scalar  moment  of 

each  of  these  dipoles  and  the  unit  vector  in  the  direction  of  the 

C-CC  tend  attached  to  the  k-th  carbon  atom.  Then 

= PaYi 

even  k 

^ — o 1 ^ 

and  /S  = )_  » 

eve* 


where 


(3‘k)=^f(bte3iy;\. 


Using  the  geometry  of  the  chain,  we  may  show  that3* 

k'  -3  “ *M*V^*i )*•••*  A(rj,Xk  A(»^, 

if  k '•*  n I . The  matrix  A(6y<p)  is  defined  by  Equation  (141).  The 
integral  for  may  be  evaluated  in  much  the  same  manner  as 

is  found  in  Appendix  II.  For  k<sn  + I , we  find  that 

U-i 

(k)_  k.o-C-'/j  (vO-,  Cl 


K -H) 


* r\-*2 


event*  q-2  »»**♦• 


2-.-I 


(229) 

(230) 

(231) 

(232) 

(233) 


*The  order  of  the  dyadics  in  Equation  (232)  is  the  reverse  of  the 
order  given  by  Kirkwood  and  Fuoss,  which  appears  to  be  erroneous 
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where  if*  is  (I  ± and  (o)  , 0)  } and  (0  are  defined  in 

t»  ^ ^ 

Appendix  II  . The  numbers  are  the  expansion  coefficients  of  stnfcJ 

in  Legendre  polynomials  of  cos  0 [17]: 


sin  e = a«  P*Cco5 ©), 

even  a 

with  ao  = TT/4-  and 


c. 


(*-  l)(«  + 0(2»  + 5)  , 

(«>2)(«+4)(2-.l)  ^ 


(234) 


(235) 


The  quantity  <x  is  equal  to  t)  s number  of  that  are  ±1: 


From  the  form 

r=^(Ov  (236) 

V2 

„oo  . , 

of  for  K ^ n + l } it  may  be  seen  that  is  zero 

unless  there 

exists  a k such  ihai 

1) 

3 

0 

11 

O 

2) 

= 0 for  q,  i n 

3) 

- «r< 1 

4) 

j rn^|  =Oot-l  n-*-2^c^4.k-l 

5) 

(0V0)v  = 0 *or  n+Ucl>  4 u'2 

6) 

|mj  = 1 

7) 

=.  0 fov  ^ ^ k + 1 . 

Corresponding  selection  rules  occur  if  k $ n , From  6)  and  7),  we  may 

-(  k) 

show  that  there  is  at  most  one  nonzero  p in  the  sum  for  jxx  , It 
follows  that 


/v;- 


1 

r ^ 


(237) 


*The  positions  of  and  if.  in  Equation  (233)  are  determined  by  the 

i-.  sign  in  Equation  (232),  Since  ^x‘k’  (3*Ch/  does  not  depend  upon  the 
choice  of  sign,  all  the  optical  isomers  of  polyvinyl  chloride  have  the 
same  dielectric  properties. 
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The  next  step  is  to  evaluate  /I CM.  The  same  approximation 
can  be  used  here  as  is  employed  in  Appendix  IIT  : 

n(A)  " r=rW  K • (238) 

N(X)  is  the  number  of  nonzero  ^ that  may  be  obtained  by  rearrange- 
ment of  the  m while  keeping  2 , k , and  <T  constant;  K is  the  average 
eigenvalue  of  these  j. ix  . For  convenience,  k will  be  replaced  by  the 
variable  s , where  s denotes  the  number  of  m between  and 

l> 

mh_,  inclusive: 

s = U-  n.  (239) 

As  an  example,  let  us  assume  that  i * I and  consequently 
that  rr>h*,=°  . From  the  theorem  in  Appendix  Til,  it  may  be  shown 

that 

MX)  = 2 (s-<r- 1 )•  (240) 

The  additional  power  of  two  appears  when  both  sides  of  the  chain  are 
included  in  W(M  . From  Equation  (209)  in  Appendix  V,  we  have 


7 --  (p  + l)(p  + 2)  D„  v - r IV"  “I 

x - 2_  I'-IVI] 

vz 

The  sum  may  be  approximated  by  an  integral,  to  give 

(p  + 2)  D,  nV  r | 


V 


A = 


2 A 


. n r r i i ' 


(241) 


(242) 


In  addition,  is  equal  to  n^/3  . Combining  Equations  (237), 


= K*  = ^L(Sj*  b.f  =n  h-^/3. 

jik  j 
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(238),  and  (240),  we  find  that  the  contribution  to  Q(X)  for  J2  = I is 

[O(A)]  = 4 H ( Z!.\  (243) 

A'i  s ,<r 

The  limits  on  this  sum  are  s s n,  s/2  i <r  < s , s even  ^ ana 


ff  4 


2XA  s 

( p ■*■  2 ) D0  n 


(244) 


If  we  introduce  the  asymptotic  expansion  of  the  binomial  coefficient, 
the  sum  may  be  replaced  by  an  integral 

im  (245) 

J 30  Ct*/i 

with  Equation  (244)  still  to  be  added  as  a limit.  The  integral  over 
<r  is  approximately  1/3  if  the  limits  include  2 s/3  , and  zero  if  they 
do  not.  Hence 


[n(AVl  = n{  *3  * 


(246) 


with 


2 _ 2*A  r i » i 

3 (p*2|0/  1(0-5)'’*'  np*'  J 


-I 


The  second  integration  gives 

[n(^]iri  = T[,“(,  + XT»)  P"l- 

As  in  Equation  (9C), 

rj  = >i2np*l/(p*2)kT. 


(247) 


(248) 


(249) 


In  order  to  find  the  total  value  of  Q(X)  , it  is  necessary 
to  add  up  the  contributions  for  each  value  of  . The  dependence  of  A 


I 
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upon  S,  can  be  neglected,  and  O(x)  - 3[n(X)]^(  : 

n(;0  = i - ( i + )~  p^1"- 


(250) 


From  Equations  (227)  and  (228)  it  is  readily  found  that 


GW  = 


.;r0  -/<>*■  (r+r;)lp*,; 


(251) 


09 

%* 


and 


0W'‘H  • 


(252) 


Ey  a change  in  variables.  Equation  (252)  becomes 


(i-vi)0(.)  - - - 

, r + i.  ujt.  - I 


(253) 


For  arbitrary  p , it  is  impossible  to  find  <3(y»)  in  closed 
form.  For  p = 5 , however,  the  integration  can  be  performed.  Ke  obtain 


0-^)0  - njj)*  (l-  J 


-I, 


(254) 


y is  defined  by 

i 

l,  * (255) 

As  above,  V±  is  . Q(uj)  is  given  in  the  table  below  for 

various  values  of  oj  , X(u>)  is  the  real  part  of  Q(.<jl>)  , while  -H(o‘) 


is  the  imaeir^ry: 


Q(o>)  = J(u>)  - L H(w). 


(256) 


! 


I 
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The  maximum  value  of 

HU) 

is  about  0.143  at 

* 10.  H (u >)} 

is  plotted  in  t igure 

6,  wit ii 

i yj  set  equal  to 

10  seconds. 

Table  III 

Theoretical  Reduced 

Polarization  for 

J(«>) 

HU) 

HU)/HW 

.01 

.997 

.0073 

.051 

.05 

.988 

.0236 

.165 

.959 

.0570 

.398 

.5 

.923 

.0837 

.584 

.7 

.904 

.0949 

.662 

X 

.880 

,1063 

.742 

2 

.825 

.1254 

.875 

5 

.740 

,1404 

.980 

7 

.707 

.1427 

.996 

10 

.673 

,1433 

1.000 

20 

.606 

.1399 

..  S76 

50 

,525 

.1290 

.900 

100 

.468 

.1193 

.632 

200 

.418 

.1085 

.756 

500 

.359 

.0947 

.661 

2,000 

,285 

.0759 

.530 

10,000 

,218 

.0583 

.407 
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